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COMPUTER-BASED MODEL FOR IDENTIFICATION AND 
CHARACTERIZATION OF NON-COMPETITIVE INHIBITORS OF 
NICOTINIC ACETYLCHOLINE RECEPTORS AND RELATED LIGAND- 
GATED ION CHANNEL RECEPTORS 

5 

The present application is a Continuation-In-Part of U.S. 
Application No. 10/411,206, filed April 11, 2003, the entire contents of 
which are hereby incorporated by reference. 

The present application includes an appended Sequence Listing of 
L0 15 amino acid sequences and Appendices 1 to 5 providing computer 
programming scripts, parameter files and atomic coordinates of computer 
models of the luminal channel portion of the ligand-gated ion channel 
subtypes. 

L 5 FIELD OF THE INVENTION 

The present invention relates to a computer system for generating 
molecular models of ligand-gated ion channels and in particular, 
molecular- models- of the inner lumen of- a ligand-gated ion channel ^and. . 
associated binding pockets. The present invention further relates to a 

20 computer system simulating interaction of the computer-based model of 
the ligand-gated channel and non- competitive inhibitor compounds for 
identification and characterization of non-competitive inhibitors and to 
inhibitor compounds so discovered. The present invention also relates to 
methods for treating various disorders related to ligand-gated ion channel 

25 receptor function. The invention also provides a way to examine 
compounds for "off-target" activity that may cause undesirable side effects 
to a desired target activity or that may represent a new therapeutic 
activity for a known compound. 

30 
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BACKGROUND OF THE INVENTION 

Ligand gated ion channels (LGICs) are currently very important f 
targets for drug discovery in the pharmaceutical industry. The 
superfamily is separated into the nicotinic receptor superfamily (muscular 
and neuronal nicotinic, GABA-A and-C, glycine and 5-HT3 receptors), the 
excitatory amino acid superfamily (glut am ate, aspartate and kainate 
receptors) and the ATP purinergic ligand gated ion channels. These 
families only differ in the number of transmembrane domains found in 
each subunit (nicotinic - 4 transmembrane domains, excitatory amino 
acid receptors - 3 transmembrane domains, ATP purinergic LGICs - 2 
transmembrane domains) . 

Nicotinic acetylcholine receptors (nAChRs) are a family of ligand 
gated ion channels that control the fast permeation of cations through the 
postsynaptic cell membrane when stimulated by acetylcholine. 
Physiologically, nAChRs are key targets in drug discovery for a number of 
diseases, including Alzheimer's and Parkinson's disease, and have been 
widely discussed and investigated. 

^Structural and functional studies of 'nAChRs have led tcf 'the 
elucidation of three physiological states of the receptor: 1) resting 
(channel closed); 2) acetylcholine stimulated (channel open); and 3) a 
desensitized state where the ion flux is inhibited even in the presence of 
neurotransmitter. The overall structure of nicotinic acetylcholine receptor 
of Torpedo marmorata has been examined by Unwin and coworkers using 
cryo-electron microscopy and revealed the conical shape of the channel 
portion of the receptor and the relationship of the membrane- spanning 
helices to each other (see Figure 1). In spite of these unprecedented 
advances in resolving the structures of transmembrane proteins, the 
detailed, atomic resolution, structure of the entire nAChR family remains 
unresolved. . 

Muscular nAChRs are located at the nerve-muscle junctions and 
are responsible for triggering motor motion, and neuronal nAChRs, widely 
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distributed in the nervous system, are involved in the fast synaptic 
transmission of inter-neuronal communication. It is known that these 
receptors are structurally similar in their overall composition but differ in 
the exact make-up of the protein subunits forming the receptors. 
5 The nicotinic acetylcholine receptor (nAChR) is presently the best 

characterized member of the ligand-gated ion channel superfamily. The 
nicotinic receptors are of great therapeutic importance. The subunits 
assemble combinatorily to form a variety of pentameric transmembrane 
protein subtypes. 

10 Each receptor is formed by bringing together five separate trans- 

membrane proteins, each containing a large extra- cellular N-terminal 
domain, four membrane spanning alpha helices (Ml, M2, M3, and M4) 
and a small C-terminal domain (see Figure 1). Two, homologous, 
neurotransmitter binding sites are formed by the N-terminal domains 

15 where cholinergic agonist and competitive antagonists bind, and are the 
usual targets for drug design. The ion channel is formed by a pentameric 
arrangement of the M2 helical segments contributed by the five proteins 
(see Figure ~2). The channel specificity, characteristic of each receptor 
subtype, is controlled by the identity of each of the M2 helices. 

20 Neuronal nicotinic acetylcholine receptors (nAChRs) are the class of 

ligand-gated ion channels of the central and peripheral nervous system 
that regulate synaptic activity. The basic structure of the nAChR is shown 
in Figures 1 and 2. Referring to Figure 1, nAChR consists of five 
transmembrane subunits 1, 2, 3, 4, 5 oriented around a central pore 6 

25 permeable to cations. Cations flow through the pore is regulated by ligand 
binding. The subunits in nAChR are typically a subunits and p subunits. 

At present, 12 different homologous subunits have been identified 
in neuronal nAChRs, 9 a subunits (a2-al0) and 3 p subunits (p2-p4). The 
major difference between a and p subunits is the presence and location of 

30 the disulfide bond formed by two adjacent cysteines in the a systems, the 
absence of this feature distinguishes non-oc subunits. This disulfide bond 
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located on the extracellular domain plays an important role in 
neurotransmitter binding as well as the mechanism of channel opening. 
These subunits combine to form multiple nAChR subtypes and 
predominant stoichiometry is (cc)2(p)3, however pentamers containing only 
5 a subunit are also known e.g., (cc7)s. In case of muscular nAChR the 
stoichiometry is more complicated, the muscular nAChR receptor is 
predominantly described as (oc^pSy. 

The nAChRs are very complex systems with dozens of potential 
different binding domains for different classes of compounds of both 

10 endo- and exogenous origin (Arias H.R., (1997) Topology of ligand binding 
sites on the nicotinic acetylcholine receptor. Brain Res. Rev. 25: 133-91). 
Two primary cholinergic binding sites are located on the extracellular side 
7 (refer to Figure 1 , approximately 30-35 A above the membrane) in the 
pocket at the interface between the a and (3 subunits. The nAChR 

15 contains several other classes of binding sites at which non-competitive 
inhibitors (NCIs) bind (Arias H.R. (1998) Binding sites for exogenous and 
endogenous non-competitive inhibitors of the nicotinic acetylcholine 
receptorr Biochira. Biophyg: Act. 1376: 173-220). One, so-called "luminal 
high affinity 5 " NCI binding domain is located on the surface of the internal 

20 lumen forming the ion channel. This site is a highly polar and negatively 
charged domain, which primarily plays the role as a cation selector. In 
general, an NCI compound does not compete with the neurotransmitter 
ligand of the receptor for binding to the neurotransmitter ligand binding 
site of the receptor located on the external surface both a subunits in a 

25 pocket approximately 30-35 A from the transmembrane portion of the 
subunit (that is, above the surface membrane when the receptor is 
expressed on in a cell), as described by Arias [Arias, H.R. (2000) 
Localization of agonist and competitive antagonist binding sites on 
nicotinic acetylcholine receptors Neurochem. Int 36, 595-645]. 

30 Such drugs as mecamylamine, ketamine, bupropion or barbiturates 

bind in the narrowest region of the channel on the cell membrane level. 
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Inhibitors acting there are mainly amines. It is believed that the ligands 
bind into this region and sterically plug the channel, blocking the flux of 
ions. 

"Non-luminaF sites are the population of 10-30 binding sites 
5 located mostly at the lipid-protein interface for which an allosteric 
mechanism of non-competitive inhibition was proposed. Agents of 
different origin (steroids, fatty acids, alcohols, local anesthetics etc.) can 
bind to those sites and modulate nAChR activity. 

Other classes of ligand-gated ion channels include GABA (Johnston 

10 G. A. (2002) Medicinal chemistry and molecular pharmacology of GABA(C) 
receptors. Curr Top Med Chem 2, 903-13), 5HT3 (D.C. Reeves, S.C. 
Lummis, (2002) The molecular basis of the structure and function of the 
5-HT3 receptor: a model ligand-gated ion channel (review). Mol Membr. 
Biol 19, 11-26), AM PA (T.B. Stensbol, U Madsen, P. Krogsgaard-Larsen, 

15 (2002) The AMPA receptor binding site: focus on agonists and competitive 
antagonists. Curr. Pharm. Des. 8, 857-72) and NMDA (K.A. Macritchie, 
A.H. Young, (2001) Emerging targets for the treatment of depressive 
disorder. Expert Qpin. Then Targets 5, 601-612) receptors, etc. Although 
the molecular structure of these receptors differ significantly, it is believed 

20 that the luminal domains are homologous to the luminal domain of 
nAChRs. There are five (or occasionally four) transmembrane helices 
forming the wall of the channel with "rings" of polar amino-acids exposed 
on the pre-forming surface and the same non-competitive inhibition 
phenomenon can be observed. 

25 In summary, the luminal high affinity NCI binding domain is 

located on the surface of the internal lumen forming the ion channel. 
Drugs of different origin bind in this region and sterically plug the 
channel blocking the flux of ions. 

Non-competitive inhibition of the nAChR can be responsible for 

30 severe adverse drug effects. On the other hand, designing ligands that 
specifically interact with this site can be part of the development of new 
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treatments of Alzheimer's and Parkinson's diseases, for example by 
identifying compounds likely to exhibit side effects through non- 
competitive inhibition of a LGIC. Furthermore, the compounds identified 
as NCIs by the present method are likely to find use in treating Tourette's 
5 syndrome and cognitive disorders, schizophrenia, pain [see, Lloyd, G.K. 
and Williams, M. (2000) J. Pharmacol Bxper. Ther. 292, 461-467.], 
anxiety, depression, neurodegeneration and addictions caused by an 
overactive LGIC receptor, especially diseases in which nicotine agonist 
activity against a neuronal nAChR is part of the etiology (e.g. smoking 

10 addiction). The invention can also be used to evaluate cardiovascular 
toxicity of a compound mediated by non-competitive inhibition of a LGIC 
receptor, e.g. arrythmia and GI spasming or diarrheal side effects of a 
compound caused by inhibition of a musclar nAChR. 

Classical methods of NCI identification are time consuming and not 

15 effective in rapid screening of chemical libraries of drug candidates. 

Several different molecular models of the nAChR transmembrane 
domain have been reported (Capener CE, Kim HJ, Arinaminpathy Y, 
Sansom MS (2002) Ion channels: structural bioinformatics and modelling. 
Hum Mol Genet 11:2425-33). However, none of those models were used to 

20 investigate interaction with channel blockers. A computer based model for 
in silico simulations of NCI interactions with the luminal domain of LGICs 
is needed to better understand the phenomenon of the receptor's 
inhibition by NCIs. 

Furthermore, in drug discovery, the potential adverse effects of drug 

25 candidates are of great importance. In-depth understanding of 
mechanistic interaction of luminal NCIs with different subtypes of LGICs, 
especially of nAChRs, is required to remove potential unwanted side 
effects at this site. In this respect, a rapid screening technology that 
would identify NCIs of LGICs, and especially of nAChRs would be greatly 

30 desired. 
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The functional determination and characterization of a NCI of a 
LGIC is very complex and time consuming. One approach is affinity 
chromatography based on immobilized receptor protein. This is a versatile 
tool for investigation of intermolecular interactions of a receptor with its 
5 ligands. The chemometric approach of affinity chromatography can be 
employed for determination of reliable relative affinities of ligands as well 
as kinetic characterization, which otherwise would be inaccessible, for a 
large set of compounds (Kaliszan R., Wainer I.W. (1997) Combination of 
Biochromatography and Chemome tries: A Potential New Research 

10 Strategy in Molecular Pharmacology and Drug Design. In 
Chromatographic Separations Based on Molecular Recognition. K. Jinno, 
editors Wiley-VCH). 

Methods using nAChR and other receptors immobilized on a 
chromatographic support have been elaborated (US patents 6,387,268, 

15 6,139,735, provisional application no. 60/337,172). It was shown that the 
obtained stationary phases worked as selective binding materials for 
competitive cholinergic ligands and can be used for high throughput 
screening of various competitive agonists and antagonists (R. Moaddel, 
LW. Wainer, (2003) Immobilized nicotinic receptor stationary phases: 

20 going with the flow in high-throughput screening and pharmacological 
studies J Pharm Biorned Anal. 30, 1715-24). The usefulness of such 
columns based on immobilized nAChR for investigations and modeling of 
NCI affinity has also been demonstrated. Using a novel non-linear 
chromatography approach off and on kinetics of ligand interaction with 

25 the receptor can be determined. (K. Jozwiak, J.Haginaka et al., (2002) 
Displacement and nonlinear chromatographic techniques in the 
investigation of interaction of noncompetitive inhibitors with an 
immobilized a3p4 nicotinic acetylcholine receptor liquid chromatographic 
stationary phase. Anal Chem 74: 4618-4624). 

30 

BRIEF DESCRIPTION OF THE DRAWINGS 
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The features of the invention may be better understood by reference 
to the drawings described below. 

Figure 1 and Figure 2 schematically show the general structure of a 
neuronal nicotinic acetylcholine receptor (nAChR). 
5 Figure 3 shows the luminal domain of a3p2 and a3p4 channels. 

Red - negatively charged (e.g., aspartic acid, glutamic acid), orange - polar 
(e.g., serine or threonine), green - hydrophobic (e.g., leucine, valine), light 
blue - positively charged (e.g., lysine), dark blue - aromatic (e.g., 
phenylalanine). Figure 3a shows a model of the a3(32 luminal domain 

10 having five helices forming the wall of the ion channel. Figure 3b shows 
the a3(32 luminal domain model (in ribbon and CPK rendering) in 
perpendicular view. Only 3 helices are shown for clarity. Figure 3c shows 
a model of the a3(34 luminal domain having five helices forming the wall 
of the ion channel. Figure 3d shows the a3p4 luminal domain model (in 

15 ribbon and CPK rendering) in perpendicular view. Again, only 3 helices 
are shown for clarity. The cleft formed by the substitution of 
phenylalanine for valine at position 15 in the helix is noted by the arrow. 

Figure 4 is a schematic representation of a computer system useful 
in the practice of the invention. 

20 Figure 5 is a model of luminal domain of a3p4 subtype of nAChR 

illustrating its electrostatic potential of the inner surface of the channel. 
The Figure particularly shows the electronegative potential of the cation 
selector region of the channel. Negative potentials are shown in red, and 
the positive potentials are shown in blue. 

25 Figure 6 shows a two cluster interaction of the ligand PCP with 

a3p4. Generally NCIs bind into the small pocket formed on the apolar 
domain (Phenylalanine /Valine rings). Tested structures primarily entered 
a hydrophobic pocket formed between the a3 and P4 helices and 
subsequently interacted with protein side chains forming hydrogen bonds. 

30 Ligands most likely form two separate clusters on two symmetrical active 
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sites. Estimated free energies of docking are in the range of experimental 
IC50 of tested inhibitors. 

Figure 7 shows example compounds tested by chromatography on 
an Q3p4 nAChR affinity column. Among the tested drugs are aliphatic 
amines like mecamylamine, amantadine, memantine and such compound 
like bupropion, ketamine and mk-801. Also, some examples of more 
complicated structures include clozapine, pep, methadone and verapamil. 
Further, the structures of two enantiomers dextromethorphan and 
levomethorphan. Finally, there is a structure for ethidium: the only 
compound permanently ionized and that binds to its specific site. 

Figure 8 shows the mecamylamine binding to the luminal domain of 

a3p4. 

Figure 9 shows the MK-801 binding to the luminal domain of a3p4. 

Figure 10 shows a correlation of log k' (chromatographic) with log 
(1/ki) (docking simulation). 

Figure 11 shows the enantio selectivity of the 
dextromethorphan/ levomethorphan pair determined in chromatographic 
experiments. Dextromethorphan had a -longer retention time and the 
profile was more asymmetric. 

Figure 12 shows an overlay of the most stable docked orientations 
of dextromethorphan (grey) and levomethorphan (magenta) complexes 
with the a3p4-nAChR luminal domain (Figure 12a) and the a3|32- 
nAChR luminal domain (Figure 12b). Only two helices (one ct3 and one P) 
are presented for clarity. The serine residues (position 8, and also in 
position 12 in Figure 12a) are rendered as CPK models and colored in 
orange. The cleft formed by the valine to phenylalanine substitution in 
the 04 subtype helix is indicated by the arrow in Figure 12a and its 
absence is similarly indicated in Figure 12b. 

Figure 13 illustrates cluster analysis according to the present 
invention. Three clusters were identified among tested NCIs. Cluster 1 - 
red, cluster, 2 - magenta, cluster 3 - green. 
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Figure 14 shows the synthetic scheme for novel compound DM-01; i 
- 1-chloroethyl chloroformate; ii - methanol; iii - 1-chloroacetone. 

Figure 15 shows the NCI activity in the Rb + efflux assay of 
compound DM-01. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention results from understanding of the 
interactions between a particular subtype of the neuronal nAChRs and 
molecules that inhibit the flow of ions through the cell membrane. A first 

10 step in this understanding is to characterize the composition of the 
membrane-spanning M2 helices. So far, twelve distinct M2 helices (known 
as subunits of the channel), nine labeled alpha (a2-al0) and three labeled 
beta (p2-(J4), have been shown to form channels of a wide variety of both 
homomeric and heteromeric subtypes of neuronal nAChRs. The most 

15 common subunit stoichiometry has been determined to be (aX)2(f3Y) 3? (X = 
2-4; Y = 2-4), respectively for heteromeric subtypes and (aZ)s ? (Z = 7-10) 
for the homomeric subtypes. However, other, more complex, combinations 
have also been reported. These various subtypes have been found in 
different locations of the central and peripheral nervous system and can 

20 be assigned to different functions. For instance: the a4(32 and a4(34 
subtypes appear to play a role in cognition, neurodegeneration, pain, 
anxiety and depression; the oe3(32 subtype in dopamine release and 
Parkinson's disease; the oc7 in GABA release; the oc9 in auditory function 
and development; and the oc3|34 in norepinephrine release, cardiovascular 

25 and gastrointestinal action. In addition, NCIs of ligand-gated ion channels 
are expected to have therapeutic benefit in treatment of cognitive 
dysfunction/ attentional disorders such as ADHD, neurodegenerative 
diseases such as Alzheimer's disease, schizophrenia, depression, epilepsy, 
Tourette's syndrome and in smoking cessation. 

30 Non-competitive inhibition of nAChRs may be responsible for many 

of the adverse effects attributed to drug therapy. For example, the 
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impairment of cardiovascular function observed during ketamine 
anesthesia has been associated with the inhibitory action of ketamine on 
ganglionic nAChRs. The administration of such drugs as methadone 
(opioid antagonist) mecamylamine or verapamil (antihypertensive agents) 
5 often results in gut motility impairment and constipation, which has been 
associated with their NCI activity on the a3p4 nAChR. 

The pentameric bundle of M2 helices forms the "lumen" (Figures 
1-3), the central surface of the narrowest part of the channel, which takes 
part in channel gating and ion selection. The amino acid residues forming 

10 the surface of the lumen are quite conserved across different subunits, 
and form distinct regions of the channel, or "rings" (see Figure 2 and 
Table 2). These rings are important for proper function and selectivity of 
the neuronal nAChRs and are common to all subtypes. An illustration of 
this importance is the fact that even a single point mutation in this 

15 domain can lead to a variety of serious diseases e.g., autosomal dominant 
nocturnal frontal lobe epilepsy, associated with a serine (S) ~> 
phenylalanine (F) mutation in the M2 segment of the oc4 subunit of 
nAChR (Steinlein, O. K. Nicotinic acetylcholine receptors and epilepsy; 
Curr. Drug Target CNS Neurol Disord. 2002, 1, 443-448.). Therefore, the 

20 sequence and structure of the M2 subunits forming the luminal domain 
are important for understanding disease states associated with nAChRs. 

The luminal domain of the ion channel has been identified as a high 
affinity binding site for a large number of exogenous and endogenous 
substances in both the open and desensitized state. Many drugs, 

25 particularly ionizable amines, can elicit deleterious side effects by binding 
to the surface of the lumen, sterically plugging the channel and blocking 
the flux of ions. This mechanism is distinct from the traditional 
cholinergic mechanism of receptor regulation, and ligands inhibiting the 
receptor in this way are called non- competitive inhibitors (NCIs) or 

30 channel blockers. Noncompetitive action on the neuronal nAChR has 
been assigned to a large number of marketed drugs and their metabolites 
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and can be responsible for many toxic side effects of various therapies. 
For example clinical side effects observed during ketamine anesthesia 
(i.e., the impairment of the cardiovascular function, etc.) have been 
associated with the inhibitory action of ketamine on ganglionic nAChRs 
5 (Friederich, P.; Dybek, A.; Urban, B. W. Stereospecific interaction of 
ketamine with nicotinic acetylcholine receptors in human sympathetic 
ganglion-like SH-SY5Y cells. Anesth. 2000, 93, 818-24.). Thus, there is a 
need to develop models to identify ligands that might be NCIs of nAChRs. 

The present invention lies in part in a computer system that 

10 generates molecular models of ligand-gated neuronal receptors and a 
method of using the same. The computer system generates a computer- 
based model of the inner lumen of a ligand-gated ion channel having 
binding pockets for non-competitive inhibitors. The computer system 
simulates interaction of structures from chemical libraries or of any 

15 desired compound with the generated computer-based model of the 
ligand-gated ion channel. The simulation can serve to predict and 
describe the pharmacological importance of the interaction. Thus, the 
-invention-constitutes a system for drug discovery and for screening of a 
drug candidate for unexpected side effects and toxicities. 

20 In an embodiment of the present invention, as shown in Figure 4, 

the computer system comprises a memory, e.g. disk 105, storing 
positional data of the atomic coordinates of the transmembrane portion of 
at least one subunit of a ligand-gated neurotransmitter receptor protein, 
and a processor 101 generating a molecular model having a three 

25 dimensional shape representative of the pore portion of the ligand-gated 
neurotransmitter receptor based on positional data. During execution of 
the process for generating the molecular model, it is understood that the 
positional data would be stored in, for example, RAM 102, or other 
memory readily accessible by the processor 101. 
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The memory, in particular, stores data of the atomic coordinates of 
at least an a chain and a p chain of a nicotinic acetylcholine receptor. The 
data of the atomic coordinates can include atomic coordinates of at least 
5 one polypeptide having an amino acid sequence selected from the group 
consisting of the polypeptides shown in Table 1 (SEQ ID NOS: 1-15). The 
data of the atomic coordinates should include atomic coordinates of the 
portion of the transmembrane portion of the subunit consisting of at least 
the amino acid sequence of residues 8 to 19 of SEQ. ID NOS: 1-15. 

10 The processor 101 can generate a molecular model of the luminal 

domain portion, especially the pore, of a ligand-gated neurotransmitter 
receptor having a subunit stoichiometiy ranging from (a)s(P)o to (a)o((3)s. 
For example, the subunit stoichiometiy can include(a)2(P)s useful for 
modeling the neuronal nACliR regulating cardiovascular and GI actions. 

15 The model may be generated with only four helices to model other LGIC 
families. 

Modeling Step: 

In~ generating a molecular model and simulating its interaction with- 
various molecules, the computer system of the present invention first 

20 generates a molecular model of the receptor channel based on a template 
structure determined in an NMR investigation of synthetic channel model 
(Opella S.J., Marassi P.M., Gesell J.J., Valente A.P., Kim Y., Oblatt-Montal 
M., Montal M., (1999) Structures of the M2 channel-lining segments from 
nicotinic acetylcholine and NMDA receptors by NMR spectroscopy. Nat 

25 Struct Biol 6:374-9). Using this model, the molecular structures of all of 
the neuronal subtypes of nAChR can be built. All subtypes of nAChR 
share several common structural arrangements in the luminal domain, 
which makes it possible to build the model of a particular subtype using a 
homology modeling approach. 

30 Once a molecular model is generated, the model is refined. A 

preferred software package for refining the molecular model is the AMBER 



WO 2004/091492 



PCT/US2004/0 10978 



[- 17 -] 

molecular modeling package, e.g. AMBER version 7, (D.A. Pearlman, D.A. 
Case, J.W. Caldwell, W.S. Ross, T.E. Cheatham III, S. De Bolt, D.M. 
Ferguson, G.L. Seibel and P. A. Kollman, (1995) AMBER, a package of 
computer programs for applying molecular mechanics, normal mode 
5 analysis, molecular dynamics and free energy calculations to simulate the 
structural and energetic properties of molecules. Comp. Phys. Comm. 91, 
1-41). The AMBER package contains a set of molecular mechanical force 
fields for the simulation of biomolecules and a package of molecular 
simulation programs. In particular, the model is preferably refined using 

10 the "SANDER" program (for Simulated Annealing with NMR-Derived 
Energy Restraints) was used. SANDER is the main program used for 
molecular dynamics simulations. SANDER allows for NMR refinement 
based on NOE-derived distance restraints, torsion angle restraints, and 
penalty functions based on chemical shifts and NOESY volumes. 

15 Once the model has been refined using the SANDER program of 

AMBER, the final model is evaluated. A preferred software package for 
evaluating the final model is the PROCHECK package, e.g. version 3.5.4 
(Laskowski R A, MacArthur.M W, Moss D S & Thornton J M, (1993) 
PROCHECK: a program to check the stereochemical quality of protein 

20 structures. J. Appl. Cryst, 26, 283-291). PROCHECK checks the 
stereochemical quality of a protein structure, producing a number of 
PostScript plots analyzing its overall and residue-by-residue geometry. 

In order to construct subtype-specific molecular models, the 
primary structures of the particular subtypes are required. Different 

>5 subtypes can be found in different region of the human brain and 
peripheral nervous system and are responsible for specific functions. 
Subtype-specific models of the nAChR luminal domain can be utilized in 
designing subtype-specific NCIs. 

The procedure of building the luminal model can be easily adopted 
;0 to constrain models of luminal domain of other subtypes of the nAChR 
and with some modification to constrain lumen models of other classes of 
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ligand-gated ion channels. The procedure is basically explained in the 
modeling step of Example 1. The model of the a3p4-nAChR can serve as 
the template to constrain other neuronal and muscular subtypes: since 
those subtypes are very homologous (see Table 1). Only a few residues 
5 need to be modified in order to obtain new subtype. The new model after 
residue modification must be subjected to energy minimization by AMBER 
procedures described previously and finally should be evaluated using 
PROCHECK. Elaborated docking procedures can be applied to those 
models and the entire approach can be used in detailed molecular 

10 characterization of the luminal domain of specific subtypes of nAChR and 
moreover, subtype specific interaction with different classes of NCIs. 

More complicated procedures must be applied if one want to obtain 
a model of the domain formed by other classes of ligand-gated ion 
channels (GABA, NMDA, 5HT3 etc). First, ammo-acid sequence alignment 

15 modeling is performed. An example and detailed description of such 
analysis can be found in the paper by Bertaccini and Trudel [E. Bertaccini 
and J.R. Trudell, (2002) Predicting the transmembrane secondary 
^structure of ligand-gated ion channels Protein Eng. 15, 443-453]. Thus, - 
homologous parts of the ion channel can be found and a new model of 

20 transmembrane domain LGIC can be made. For some LGICs, the 
transmembrane domain is formed by four transmembrane helices instead 
of five as in the case of nAChR. In such case one of the helices must be 
removed and the remaining four need to be properly repositioned in order 
to form the channel structure. Then the model can be relaxed and refined 

25 in AMBER procedures and finally evaluated in PROCHECK. In case of 
such distinct models the docking procedures need be parameterized by 
initial studies as described in the simulation step of Example 1 below. 
The values of the size of the grid box, the dielectric constant and the 
ga_num_evals must be optimized, since the size and environment of the 

50 channel would have been changed significantly. 
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Using the modeling method of the invention, it has been discovered 
that there are NCI binding sites at the interface between a and p helices of 
LGICs, especially of nicotinic AchRs. Among modeled candidate NCIs, the 
compound enters into a small hydrophobic pocket formed by residues 12, 
5 15 and 18 of the transmembrane domains of the receptor subunits (e.g. 
SEQ ID NOS: 1-15, Table 1). A hydrophobic group of the NCI compound 
will interact with this portion of the NCI binding site. A polar group (e.g. 
an amino group) of a putative NCI can interact by hydrogen bonding with 
surrounding polar residues (e.g. residues 12 andl4 of SEQ ID NOS: 1-10). 

10 Simulation step: 

After generating the molecular model, the final molecular model is 
used as a target protein for docking simulation for compounds that may 
be potential inhibitors. A preferred software package for docking 
simulation is the AutoDock package, e.g. version 3.5. AutoDock allows 

15 docking of a flexible ligand into a rigid structure of the target protein 
using genetic algorithms as the search method. 

A particular genetic algorithm included in the AutoDock package is 
the^Lamarckian genetic algorithm. The Lamarckian genetic algorithm was — - 
preferably used with local search in order to improve efficiency. The 

20 Lamarckian genetic algorithm works in a reverse order compared to 
typical genetic algorithms. In particular, new traits in an organism 
develop because of a need created by the environment and these acquired 
characteristics are transmitted to its offspring. In AutoDock the ligand's 
atomic coordinates represent a genotype and fitness is represented by 

25 interaction free energy with the proteins. Genotypes are found through 
interations of the local search and then the atomic coordinates are 
translated into the ligand's state coordinates as the phenotype. In other 
words, in AutoDock local search is used to update the fitness associated 
with an individual in the genetic algorithm selection. 

30 The Lamarckian genetic algorithm uses as input a grid data set 

produced by the AutoGrid module. The AutoGrid module is used to create 
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3-dimensional maps over the host protein using several atom specific and 
electronic probes at each grid point. 

Results of these simulations allow the classification of tested 
compounds in terms of free energy of binding, which leads to the 
5 identification of ligands that may be potent inhibitors. The same approach 
can be used to design new compounds with high affinity binding 
properties to a specific subtype of the nAChR. A compound that is 
identified as a non-competitive inhibitor of a LGIC is one having a AG less 
than -6 kcal/mol, preferably less than -7 kcal/mol, still more preferably 

10 one having a AG less than -10 kcal/mol. 

The ligand structures used in docking simulations are preferably 
made using the HyperChem package (of HyperCube Inc., Gainsville, FL). 
In particular, it is preferred that the AMI semiempirical method 
implemented in HyperChem be used to minimize the system energy and 

15 to calculate atomic charges in final structures (J. J. P. Stewart, 
Semiempirical molecular orbital methods, in: K.B. Lipkowitz, D.B. Boyd 
(Eds.), Reviews in Computational Chemistry, vol. 1, VCH, New York, 1990, 
pp. 45=81) 

The in silico approach described above can be supported by 
-0 examining the NCI-nAChR interaction by affinity chromatography 
(Jozwiak K, Haginaka J, Moaddel R and Wainer IW (2002) Displacement 
and nonlinear chromatographic techniques in the investigation of 
interaction of noncompetitive inhibitors with an immobilized nicotinic 
acetylcholine receptor liquid chromatographic stationary phase. Anal 
55 Chem 74: 4618-4624), preferably in an iterative fashion. Chromatographic 
affinity screening can provide experimental data that is then employed for 
proper parameterization of the computer-based molecular simulation. 
Alternatively, the results of computer-based simulation can be related and 
evaluated by further chromatographic and functional experiments. 
10 Until recently, the screening of drug candidates for activity as NCIs 

was not a standard procedure in the drug development process. However, 
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the present invention will permit pharmaceutical companies to rapidly 
screen their potential drugs for NCI properties. In addition, the luminal 
domain of nAChR can be used as a target in drug discovery programs, 
which represents a new therapeutic approach to the treatment of diseases 
5 such as Alzheimer's and Parkinson's diseases and for treatment of drug 
and tobacco dependency, which are related to LGIC functions, especially 
to nAChR functions. 

The nAChR, for example, was found to contain two cholinergic 
agonist binding sites located at the interface between the a and 0 

10 subunits and on the extracellular N-terminal of the a subunits. These 
sites are key targets for drug discovery in a variety of diseases, including 
Alzheimer's disease (a 4 p2)> Parkinson's disease (a.302), cardiovascular and 
GI actions (a 3 p 4 ) ? anxiety and depression (a 4 (34), short term memory (a7) 
and auditory function and development (a9). 

15 Candidate NCI compounds discovered by the computational 

modeling method of the invention can be confirmed by in vitro 
experimental methods. Two preferred methods are by binding experiments 
or by functional assays. Either of these methods may employ the target- 
LGIC, a population of LGICs representing the target receptor and 

20 receptors that the compound should preferably not inhibit (to avoid side 
effects), or a population of LGICs representing a group of target receptors 
(with or without a group representing LGICs that the compound should 
preferably not inhibit). The LGICs for the in vitro functional assays can be 
present either as expression products in cells, as partially purified 

25 proteins, e.g. membrane preparations made as known in the art, or as 
isolated proteins. If isolated proteins are used in binding experiments, the 
proteins are preferably immobilized. 

A preferred binding assay is a displacement assay performed as 
described by Jozwiak et al. [Jozwiak K, Haginaka J, Moaddel R and 

50 Wainer IW (2002) Displacement and Nonlinear chromatographic 
techniques in the investigation of interaction of noncompetitive inhibitors 
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with an immobilized oc3p4 nicotinic acetylcholine receptor liquid 
chromatographic stationaiy phase. Anal Chem 74: 4618-4624.] Using this 
assay, a compound is identified as a non- competitive inhibitor of the 
ligand-gated neurotransmitter receptor as one that specifically binds to 
5 the ligand-gated neurotransmitter receptor with a k ? value greater than 8, 
preferably with a k ? value greater than 9 or even more preferably a k' 
value greater than 10. 

Specificity of NCI binding to particular LGICs can be shown by 
displacement of compounds that are selective for the pore portion of the 

10 desired LGIC. Specificity of the binding to a nicotinic AChR and 
homologous receptors can be shown by displacement by mecamylamine. 
Displacement of mecamylamine at a concentration of 10 \jlM (of 
mecamylamine) indicates good specific binding, ability to displace 
mecamylamine at a concentration of 40 \dM indicates strong specific 

15 binding. Preferably it is possible to displace mecamylamine at a 
concentration of 100 jiM. Thus, a compound that is a preferred NCI of a 
nicotinic AChR is one that exhibits a k' of greater than 8 in a 
chromatographic binding experiment and can be displaced by - 
mecamylamine at a concentration of 10 to 100 juIVL 

20 Preferred functional ion channel activity assays are described by 

Hernandez et al. [Hernandez SC, Bertolino M, Xiao Y, Pringle KE, Caruso 
FS and Kellar KJ (2000) Dextromethorphan and its metabolite 
dextrorphan block cc3p4 neuronal nicotinic receptors. J Pharmacol Exp 
Ther 293: 962-967] and by Jozwiak et al. [K. Jozwiak, SC Hernandez, KJ 

25 Kellar, IW Wainer (2003) The Enantioselective Interactions of 
Dextromethorphan and Levomethorphan with the a3p4-Nicotinic 
Acetylcholine Receptor: Comparison of Chromatographic and Functional 
Data submitted to J Pharmacol Exp Ther]. In brief, 1-ml aliquots of cells 
in growth medium were plated onto 24-well plates coated with poly(D- 

30 lysine). The plated cells were grown at 37°C for 16 to 18 h until reaching 
90 to 100% confluence. On the day of the experiment, the growth 
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medium was aspirated and the cells were incubated in fresh medium 
containing 2 pCi/ml 86 RbCl for 4 h at 37°C. After this loading procedure, 
the medium was aspirated and the cells were washed three times with 1- 
ml aliquots of buffer (15 mH HEPES, 140 mM NaCl, 2 mM KC1, 1 mM 
5 MgS0 4 , 1.8 mM CaCh, and 11 mM glucose at pH 7.4) to remove 86 Rb + free 
in the medium. After washing, 1 ml of buffer with or without the drugs 
under study was added to each well, and the cells were incubated for 2 
min. The incubation buffer was then collected, after which the cells were 
lysed in 0.1 N NaOH. The radioactivity in the buffer samples and cell 

10 lysates was measured by liquid scintillation counting. The total 86 Rb + 
loaded into the cells (after washing) was calculated as the sum of the 
buffer samples and the cell lysates from each well, and the amount of 
86 Rb + efflux was then expressed as a percentage of the total 86 Rb + loaded 
(fractional release). Stimulated efflux was defined as the difference 

15 between efflux in the presence and absence of nicotine (i.e., total minus 
basal efflux). The maximum 86 Rb + efflux, found at a nicotine 
concentration of -300 pM or higher, was -45% of the amount loaded and 
was independent of the amount of 86 Rb + loaded- into the cell. In studies to 
determine the inhibition of nicotine-stimulated 86 Rb + efflux by the drugs 

20 under study, data were expressed as a percentage of control values 
measured with 100 nicotine. 

In these assays a compound is identified as a NCI that inhibits the 
ion channel activity of the ligand-gated neurotransmitter receptor in 
nicotine stimulated 86 Rb + efflux with an IC50 lower than 50 jjJVL A more 

25 preferred NCI compound is one that inhibits ion efflux with an IC50 lower 
than 5 juM. Even more preferable compounds are those that inhibit ion 
efflux with an IC50 lower than 500 nM. One of skill in the art will 
recognize that compounds that are effective at even lower concentrations 
are still more preferable, and IC50 of 50 nM, or even 5 nM might be 

30 observed. 
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In some instances as described above, it might be preferred to have 
a NCI that is selective for a particular LGIC. By "selective" is meant that 
the NCI inhibits the target LGIC with an IC50 that is at least 5-fold higher 
than the IC50 of the one or more LGICs that it is desired not to inhibit. The 
5 degree of selectivity is preferably 10-fold, more preferably 20- to 50-fold, 
and still more preferably 100- to 500-fold or more. 

On the other hand, the binding assays or functional assays also can 
be used to provide initial data that can be used to constrain the in silico 
modeling method described above. Alternatively, the in silico modeling and 

10 the in vitro assays can be run iteratively to converge upon NCI compounds 
that have desired properties. 

Methods for synthesis of compounds of the invention are considered 
within the skill of the ordinary synthetic chemist. Preferred NCI 
compounds have the above structural features and exhibit activity of 

15 inhibiting the ion-channel activity of a ligand-gated neurotransmitter 
receptor in nicotine stimulated 86 Rb + efflux with an IC50 lower than 100 
jaM or other activities as set forth in detail above. 

-Dosage of compounds used for treatment of a subject can be easily 
determined by the ordinarily- skilled pharmacologist using known 

20 pharmacokinetic and pharmacodynamic assays and calculations from 
IC50 data obtained by the inventive method. Doses of from 100 |ng to 500 
mg per dose are typical. Formulation and administration of compounds 
useful for treatment is also well-known in the art. For example, many of 
the compounds listed in Table 2 have been administered therapeutically 

25 and it is expected that compounds of the invention can be similarly 
formulated and administered. 
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Example 1: Modeling of the Lumen of a q3B4 nAChR and docking of a 
putative NCI 

The molecular model of a S-M2 - nAChR transmembrane channel 
determined by frozen state NMR was used as the template for further 
modification (atomic coordinates were found in Protein Data Bank - PDB 
id: 1EQ8). This model represents a channel that mimics the 
transmembrane arrangement of known LGICs (Opella S.J., Marassi F.M., 
Gesell J.J., Valente A.P., Kim Y., Oblatt-Montal M., Montal M., (1999) 
Structures of the M2 channel-lining segments from nicotinic acetylcholine 
and NMDA receptors by NMR spectroscopy. Nat Struct Biol 6:374-9). The 
model channel consisted of 5 uniform polypeptides oriented around a 
central pore. The amino-acid sequence of this polypeptide is analogous to 
the sequence of transmembrane M2 segment of 5 subunit of nAChR found 
in Torpedo californica. 

In the 8-M2 - nAChR transmembrane channel, the spatial 
arrangement of polypeptide helices conserves five-fold symmetry, with 
certain residues exposed to the center of the pore. These residues 
(predominantly polar) form sx\ explicit surface of the channel. This is 
consistent with the concept of the presence of amino acid rings 
distributed along the pore and is a common property found in all 
subtypes of nAChR and also other ligand-gated ion channels [Changeux 
J.P., Galzi J.L. ? Devillers-Thiery A., Bertrand D. ? (1992) The functional 
architecture of the acetylcholine nicotinic receptor explored by affinity 
labelling and site-directed mutagenesis. Q. Rev. Biophys. 25: 395-432]. 

With respect to the spatial arrangement of five helices in the 
luminal domain, distribution of certain amino-acid rings along the 
channel is a common property of all subtypes of nAChR. Since primary 
sequences across different subtypes are predominantly homologous as 
presented in Table 1, and essential (exposed) residues are highly 
conserved, a subtype specific model of the luminal domain can be built 
using homology modeling techniques. 
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Based on the sequence comparison presented in Table 1, the initial 
model was modified by exchange of 5 helix residues into a3 and (54 using 
the SYBYL 6.8 molecular modeling system (Tripos Inc., 1699 South 
Hanley Road, St. Louis, Missouri, 63144, USA). Therefore, the channel 
5 containing oc3, (34, cc3, p4 and p4 helices, respectively, was constrained. 

The model was further refined by energy minimization using the 
Sander_Classic module of AMBER 6.0 software. Both termini of each helix 
were blocked in a standard AMBER procedure: acetyl beginning groups 
(ACE) and N-methylamine ending group (NME) groups were attached, 

10 respectively, to each helix. The AMBER '94 force field (Cornell, W.D., 
Cieplak, P., Bayly, C.I., Gould, I.R., Merz, Jr. K.M., Ferguson, D.M., 
Spellmeyer, D.C., Fox, T., Caldwell, J.W., Kollman, P.A., (1995) J. Am. 
Chem. Soc. 117, 5179-5197) parameters were used for energy 
minimization with the convergence criterion of the root-mean- square of 

15 the gradient to be less than 1.0E-4 kcal/mole A. Each minimization run 
was started with the steepest descent followed by the conjugate gradient 
method. A distance-dependent dielectric function was used to evaluate 
the electrostatic energy. The energy minimization run was carried-out in 
stages by relaxing i) only hydrogen atoms, ii) hydrogen + side-chain 

20 atoms, or iii) all atoms except alpha-carbons. Finally, a restrained 
minimization was also performed on the alpha-carbons of all the 
chains /residues of the model. This was to relax the structure but keep it 
near the initial position of the known template structure (PDB accession 
no. 1EQ8). Respective scripts used to run model refining with AMBER are 

25 presented in Apendix 1. 

Using PROCHECK to evaluate the model it was found that the 
whole luminal domain is constrained fully by a-helix secondary structure. 
Along the lumen model seven rings of residues exposed to the center of 
the channel can be found; three polar residues (E, T and S) and then 

30 three apolar residues (L, V/F and LL) and the last polar residue (E/K). 
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It is believed that apolar rings in the middle of the structure form 
the actual "gate" of the channel and play a role in conformational change 
of the receptor from a closed to an open state. Polar residues on both 
sides of the "gate" participate in the cation selective function of the 
5 receptor. An important structural parameter found in the obtained model 
is the change in position from valine in the a3 sequence to phenylalanine 
in the 04 sequence (see residue 15 in Table 1). This provides the formation 
of small pockets between cc3 and p4 subunits, found during the 
simulation of NCI-a3p4-nAChR interactions. The developed model of 

10 a3p4-nAChR luminal domain can be used as a template to constrain 
homologous systems of other nicotinic receptors, especially neuronal 
nicotinic receptors, and other ligand-gated ion channels. 

The resulting atomic coordinates represent the final model. Figure 5 
illustrates the electrostatic potential of the inner surface of the ion 

15 channel, and especially the electronegative potential of the cation selector 
of the channel. Figure 3c shows an example of the resulting luminal 
domain model having five helices forming the wall of the ion channel. 
Figure 3d shows a luminal domain model in perpendicular view""" with 
residue rings. 

20 In order to perform docking simulations, the Auto Grid module was 

first used to create 3-dimensional maps over the host protein using 
several atom specific and electronic probes at each grid point. An example 
parameterization file for the AutoGrid module used in this example can be 
found in Appendix 2. The optimal size of constrained grid maps was a 

25 22.5x22.5x45 A box (i.e., a grid of 60x60x120 points, each separated by 
0.375 A). This allowed exploration of the whole internal space of the 
lumen domain but prevented ligands from being bound on the external 
side. The grid-box size can bp altered in the 3 rd dimension (along the 
lumen) in order to explore interaction with a particular segment of the 

30 lumen or to calculate the interaction profile along the model. 
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An important parameter to properly explore electronic interaction in 
ligand receptor complexes is the dielectric constant value [d) used to 
calculate the electronic grid map. During the initial evaluation tests, the 
standard distant-dependent dielectric constant did not produce proper 
5 results: the electrostatic interaction were almost zero. The simulation did 
not discriminate between neutral and protonated ligands. 



Table 2. AG (kcal/mol) of best docked conformation obtained in different 
dielectric environment. 
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10 A detailed test of several d values was carried out using three pairs 

of ligands and the results are presented in Table 2. Table 2 shows an 
unexpected diminished difference between neutral and protonated 
systems when distant-dependent d was used; differences gradually 
increase with decreasing d. Simultaneously the increase in electrostatic 

15 impact in the ligand receptor interaction was noticed when a low dielectric 
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value was used. However, a very low value (d < 10) produced unrealistic 
AG values. Finally, as a mater of compromising these two effects, d = 15 
was chosen for final calculations as the value producing suitable 
electronic properties of the ligand-receptor complex in the transmembrane 
5 ion channel system. This approach is in agreement with values of the 
dielectric constant in transmembrane pores obtained by theoretical 
calculations (Cheng et al., (1998) Eur. Biophys J., 27 105-112 and 
Gutman et al., (1992) Biochim. Biophys Acta 1109: 141-148) where it was 
found that the actual dielectric constant in transmembrane channels 
10 remains low and ranges from 25 to 5 depending on the structure. Thus, in 
the case of the NCI-nAChR docking simulations d value can vary from 10 
to 20. 

The resulting ligand 3D structure was loaded into the AutoDock 
system and was iteratively sampled over previously created grid-maps in 
15 order to find optimal positions and the lowest energy of interaction. An 
example parameterization file for the AutoDock module used in this 
example can be found in Appendix 3. 

The Lamarckian genetic algorithm with local search was used from- - — 
the AutoDock package. Atomic coordinate files of ligands were 
20 transformed into a format suitable to AutoDock using the HIN2PDBQ 
script (Johansson M. (2002) Some computational chemistry related 
python conversion scripts. See Web site 

helsinki . fi / %7Empj ohans / python / ) . 

The ligand structures used in the docking simulations were made 
25 using the Hyper Chem software package. Further, the AMI semiempirical 
method implemented in HyperChem was used to minimize the system 
energy and to calculate atomic charges in final structures. 

An initial simulation was performed in order to optimize the docking 
settings. Since previously described docking space seemed to be large in 
30 the model of a304-nAChR active site (22,781.25 A 3 ) it was important to 
optimize the maximum number of energy evaluations (ga_num_evals) 
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required in each search run. It was found that too low a value of 
ga_num_evals could result in finishing the simulation too quickly, and the 
global minimum of the complex conformation may not be found. A set of 
test simulations on several ligands including conformationally flexible and 
5 rigid systems was performed. It was found that a ga_num_evals value of 
at least 5 million is required to assure obtaining a statistically significant 
number of lowest energy complexes. In the case of bigger ligand molecules 
with more than 2 rotatable bonds, the optimal value should be at least 50 
million. Higher values are acceptable; however higher values may 

10 dramatically increase the time of each simulation. 

The optimal number of docking search runs was found to be 50. 
Again the number of docking search runs can be higher, but would take 
more time for simulation and have no effect on the final result. 

AutoDock 3.5 implemented a free-energy scoring function that is 

15 based on a linear regression analysis, the AMBER force field, and a large 
set of diverse protein-ligand complexes with known inhibition constants 
(e.g. see Web site at scripps.edu/pub/olson-web/doc/autodock/). This 
function was employed to estimate the free energy change of the NCI- 
nAChR complex and eventually lead to an estimated inhibition constant of 

20 a particular ligand. Docking simulations allow quantitative classification 
of the stability of the NCI-nAChR complexes formed by tested ligands in 
terms of free energy of binding, which eventually lead to the identification 
of ligands exerting potent inhibitory properties. It was found that 
molecular systems forming the complex with AG value lower than -6.0 

25 kcal/mol should be considered as potential NCIs. Lower AG values 
represent more potent NCI compounds. Preferred NCI compounds exhibit 
a AG value lower than -7.0 kcal/mol; more preferred compounds exhibit a 
AG value lower than -10.0 kcal/mol. 

Detailed exploration of the spatial arrangement of ligand-receptor 

30 conformations leads to building a pharmacophore model of a subtype 
specific NCI-nAChR. Simulations on the a3p4 model showed that NCIs 
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bind predominantly in the channel in the apolar domain (F/V ring). 
Tested structures primarily entered a small hydrophobic pocket formed 
between ot3 and p4 subunits and subsequently interacted with protein 
side chains, forming hydrogen bonds. It is expected that this is a type of 
5 interaction that would not be found in those receptor subtypes that lack 
the bulky phenylalanine residue in this position. Since there are two 
quasi- symmetrical pockets between a3 and P4 helices in the model, 
ligands most likely form two separate clusters on these two symmetrical 
sites (e.g, Figure 6) at which the energy of interaction does not 
10 significantly differ. Estimated free energies of docking are in the range of 
experimental IC50 of tested inhibitors and also can be related to 
experimental affinity chromatography results. The model can be applied 
to a variety of compounds and is useful for in silico designing of new 
drugs with particularly high non-competitive inhibitory activity. 

15 

Example 2: Chromatographic assay of NCI activity 

Chromatographic studies based on immobilized nAChRs were 

- performed- to characterize^ ligand binding for broad groups of compounds: 

In order to further understand the mechanistic action of NCIs on the 

20 molecular level, the model of the transmembrane domain of the a3p4 
nAChR was built and used for computer simulations of docking inhibitors 
into the receptor. The entire approach allowed the classification of NCIs in 
terms of their functional effectiveness. 

Figure 7 presents compounds tested on an oc3p4 nAChR column. 

25 The chemicals can be divided into several subgroups. The first group 
contains drugs from different origin, which are well known as non- 
competitive inhibitors of nAChRs. The second group is of the 
dextromethorphan family, levomethorphan, dextromethorphan and its 
analogues, and the final group is verapamil, its congeners, and 

30 metabolites. In order to properly assess the influence of non-specific 
retention, five other chemicals (acetanilide, acetaminophen, 2,4- 
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dinitrobenzoic acid, 3,4-dimethoxybenzoic acid and phenylbutazone) were 
tested as negative controls. The affinity of ligands was investigated by 
non-linear chromatography on an a3p4 nicotinic receptor affinity column. 
106 Cells from the KXcc3p4R2 cell line were suspended in Tris-HCl 
5 [50 mM, pH 7.4] (buffer A), homogenized for 30 sec, and centrifuged at 
35 ? 000 x g for 10 min at 4 °C. The pellet was resuspended in 2% cholate 
in buffer A and stirred for 2 h. The mixture was centrifuged at 35,000 x g 
for 30 min, and the supernatant containing a3p4 nAChR-cholate solution 
was collected. 200 mg of the I AM stationary phase was added to the a3p4 

10 nAChR-cholate solution. Subsequently the solution was stirred for 1 h. 
The suspension was dialyzed against 2 x 1L buffer A for 24 h at 4 °C. The 
IAM liquid chromatographic support containing the a3(34-nAChR was 
packed into a HR5/2 glass column to form a chromatographic bed of 20 
mm x 5 mm i.d. The a3(34-nAChR column was then placed in the 

15 chromatographic system and used. 

Aqueous solutions [10 [iM] of each compound were prepared and 20 
|iil aliquots were injected into column. The mobile phase was composed of 
ammonium acetate [10 mM. pH 7.4] modified with methanol in the ratio 
85:15 (v/v). The flow rate was 0.2 ml/min and the experiments were 

20 carried out at ambient temperature. 

DM and LM were monitored in the positive ion mode (ESI+). The 
compounds were detected using single ion monitoring at m/z = 272 
{[MW+H] + ion}. The chromatograms were recorded and processed using 
MassLynx v. 3.5. (Micromass) . 

25 The non-linear chromatography approach was used to determine 

kinetics of the NCI-nAChR interaction in affinity chromatography studies. 

The mathematical model assumes limited (and a relatively low) 
number of active sites on the column. Slow association and dissociation of 
the drug-protein complex are the main cause of band broadening and 

30 asymmetry of the peak profile. The chromatographic peak profiles were 
analyzed using PeakFit v4.1 1 for Windows Software (SPSS Inc., Chicago, 
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IL). The mathematical approach used was the non-linear chromatography 
(NLC) model derived from Impulse Input Solution [Wade J L, Bergold AF 
and Carr PW (1987) Theoretical description of nonlinear chromatography, 
with applications to psychochemical measurements in affinity 
chromatography and implications for preparative-scale separations. Anal 
Chem 59:1286-1295.] and described by Equation 1 (PeakFit User's 
Manual, p. 8-25): 



1 — exp 



a, 



a 



2 J 



j ' 



a 0 



exp 



a-, 





C \ 




( „ > 




1-T 


flj X 










1 — exp 








K a 2 ' a 2) 









where: 

y - intensity of signal, 

x - reduced retention time, 

TO, v) = exp(-v) |exp(-0I 0 {l4vt)dt 



Eqn. 1 



Io( ) and Ii( ) are Modified Bessel functions 



area parameter, 



a x - center parameter, reveal to true thermodynamic capacity factor, 

a 2 - width parameter, 

a 3 - distortion parameter. 

Experimental chromatograms obtained by single injection of ligand 
into the chromatographic column with immobilized receptor were 
processed with PeakFit v4.11 software. After standard linear baseline 
subtraction, each peak profile was fitted to the NLC function. The set of 
NLC parameters (ao, ai, a2 and as) was collected for each profile and used 
for the calculation of descriptors of the kinetic interactions with the 
immobilized nAChR, dissociation rate constant (koff); equilibrium constant 
(K a ); association rate constant (Abn) real thermodynamic capacity factor 
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(&'), according to the following equations: 
k'= a x Eqn. 2 




Eqn. 3 




Eqn. 4 



Eqn. 5 



where: t 0 is the dead time of a column (time needed by non-retained 

substance to reach the detector); C 0 is a concentration of solute injected 

multiplied by a width of the injection pulse (as a fraction of column dead 
volume) . 

Thus, by analyzing the ligand in an immobilized receptor system 
four descriptors can be collected: retention (k% association rate constant 
(k on ), dissociation rate constant (k D ff) and equilibrium constant (logK). It 
was found that ligands which are non-competitive inhibitors have k' 
greater than 8, k on greater than 10 x 10" 6 M^s- 1 (preferred inhibitors have 
kon of greater than 15 x 10-6 IVHs-i especially potent inhibitors have kon 
greater than 30 x 10-6 M-is- 1 ), k ofr smaller than 15 S" 1 (preferably lower 
than 2 S" 1 ) and logK greater than 5.9 (preferably greater than 6.5). 

The kon value obtained in chromatographic experiments is the one 
which is closely correlated with IC50 values from functional in vitro or in 
vivo experiments. In the docking simulation, it is preferred that AG be 
lower than -6 kcal/mol (preferably less than -7 kcal/mol, most preferably 
less than -10 kcal/mol). In functional nicotine stimulated Rb+ efflux 
experiments, the IC50 value is preferrably lower than 100 |j.M (preferred 
inhibitors exhibit an IC50 < 10 jaM). 



Table 3. Detailed chromatographic characterization of tested non- 
competitive inhibitors k' - retention capacity factor, kon and k off are 
association and dissociation constant rates, respectively (kinetics of 
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formation and disformation of the complex in chromatographic system), 
logK is chromatographic equilibrium constant. 
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Values of log^T and k 7 presented in Table 3 can be regarded as a 
measure of relative affinity of tested NCI compounds for the nicotinic 
AChR. Among tested compounds, ethidium, clozapine, verapamil and 
5 some of its congeners (PR-22, nor- verapamil and galapamil) have the 
highest affinities towards the ct3p4 nicotinic receptor column as reflected 
by both log AT and k\ Both verapamil and nor- verapamil were tested for 
enantioselectivity of binding towards nicotinic affinity column but 
chromatographic experiments as well as NLC data did not exhibit 
10 noticeable differences between enantiomers. Interestingly, 
dextromethorphan exhibited markedly increased affinity compared to the 
optical enantiomer levomethorphan. 

The NLC approach allows estimating the kinetic rates of the 
complex formation and dissociation, fen and 7cbff, respectively. The well- 
15 known and potent NCIs mecamylamine, ketamine, ethidium and 
bupropion had high association constant rates. Ketamine, 
methamphetamine, amantadine and mecamylamine dissociated markedly 
-quicker th&rx other tested ligands. The lowest "dissociation constant rates " 
were exhibited by ethidium, clozapine and verapamil congeners. 

20 

Table 4. QSAR models build on of chromatographic data 
In general, it was found that different molecular properties correlated with 
chromatographic experimental characteristics: non specific bulkiness 
(Vol) and lipophilicity (logP; TASA) descriptors, specific shape descriptors 

25 (like Xiength), another group were electronic properties: ability to form 
hydrogen bonds, Ehomo and iVorder. Two latter descriptors can be 
associated with ability of a ligand molecule to protonation. Both 
bulkiness/ shape and H bond /protonation seem to be important in respect 
of known mechanism of inhibition by luminal NCIs: ligand must enter the 

30 polar pore, interact with negative and polar surface (primarily designated 
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for cation selection) and eventually block or inhibit the flux of ion during 
receptor's open stage. 
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Examples of complexes resulting from simulations are provided in 
Figures 6, 8 and 9. Figure 6 shows a two cluster interaction of the ligand 
PCP with a3p4. Figure 8 shows the mecamylamine luminal domain of 
5 a3p4. Figure 9 shows the MK-801 luminal domain of oc3p4. 

Quantitative results of simulated docking affinities were related to 
experimental results from chromatographic studies. Using AutoDock's 
scoring function, estimated inhibition constant were calculated. These 
values exhibited very good correlations with affinity data from NLC 
10 calculations (Figure 10). This correlation can be illustrated by equation: 

log *'= 0.41 8(±0.037) log(l IK.)- 0.89(±0. 1 9) 
r- 0.930 F = 127.7 n = 22 

Table 5. the collection of dextromethorphan (DM) / (LM) levomethorphan 
15 characterization by different approaches (chromatographic and docking 
were explained above), functional in vivo is nicotine stimulated Rb+ efflux 

.experiments: it was found that DM has significantly longer recovery time- 

than LM, which was predicted by chromatographic and docking modeling. 
The IC50 does not significantly differ. 
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Enantiomers have identical physiochemical properties and, therefore, all 
5 possible non-specific interactions between the enantiomers of a chiral NCI 
and an immobilized nAChR stationary phase should be equivalent. Any 
differences in the chromatographic retention between the enantiomers will 
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be due to specific binding interactions with the active site of the protein. 
Figure 11 shows chromatograms of dextromethorphan (DM) and its 
enantiomer - levomethorphan (LM). The pair of enantiomers was further 
investigated by chromatographic, docking and functional studies (Table 
5 5). It was learned from the chromatographic experiments that the drug 
dextromethorphan (DM) exert higher affinity on oc3p4-nAChR than its 
enantiomer levomethorphan (LM) and the difference in AG of the 
complexes was 0.3 kcal/mol. These data were valuable in evaluating 
parameter selection during initial tests of the docking simulations to 

10 optimally choose the channel dielectric constant or evaluate the 
usefulness of the scoring function for calculating estimated AG 
implemented in AutoDock. The docking simulations give insights into 
chiral recognition on the molecular level (Figure 12). In binding to the 
a3p4 luminal domain, both molecules interact initially with a hydrophobic 

15 pocket on the border between the a3 and (34 helices (Figures 12a). This 
binding determines the positions of the terminal amine group (blue) 
differently for dextromethorphan (grey) than levomethorphan (magenta). 
The amine group of dextromethorphan can easily form secondary 
interaction hydrogen bonds with neighboring polar residues (orange 

20 balls), while levomethorphan is less likely to form such interaction. This 
makes a difference in stabilities of two complexes by ca. 0.3 kcal/mol 
determined by both docking and chromatographic analysis (Figure 11). 

Furthermore, the estimated inhibition constant obtained during the 
simulations is very well correlated with equilibrium measures obtained in 

25 affinity chromatographic experiments. 

Example 3: OSAR - 3D clustering technique 

Classical methods for the identification and characterization of non- 
competitive inhibitors to ligand gated ion channels are time consuming. 
30 They are not applicable to the rapid screening of chemical libraries for 
potential new drug candidates nor can they be routinely used in the new 
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drug development process. An important advancement in the area is the 
development of a method of identification of potent NCIs. The method is 
based of the chemometric processing of the chromatographic data 
obtained using a stationary phase modified by immobilization of 
particular subtype of the receptor. The non-linear chromatography 
approach allows description of the NCI-receptor interactions in terms of 
real thermodynamic capacity factor (/c'), equilibrium constant for binding 
{Ka) and kinetics rate constants for association (fca) and dissociation (/ca). 
We have determined that a strong correlation exists between the drug kd 
parameter obtained in affinity chromatography experiments and the 
relative length of the effect of this drug in functional studies (nicotine 
stimulated efflux of 86 Rb + , from cells expressing the target nAChR) (K. 
Jozwiak, J. Haginaka et al, Anal Chem., 2002, 74, 4618-4624. and K. 
Jozwiak, S.C. Hernandez et al., J. Chromatogr. B. 2003, 797,423-431). 

A strong relationship between the chromatographic rate constant 
and the length of the functional effect was found. However, more than 
chromatographic affinity has been found necessary to predict the IC50 
value for NCI activity. The non-linear chromatographic parameters - 
determined in these studies were obtained in a dynamic system but under 
simplified conditions when compared to a functional assay (i.e. no 
neurotransmitter stimulation, no transmembrane potential, etc.). Thus, 
the efficacy of the NCI's expressed as IC50 values were not directly 
correlated with the calculated non-linear chromatographic parameters. 

Quantitative Structure- Activity Relationship (Q S AR) analysis 
provided models of the chromatographic affinity (Table 6). Each of the 
derived equations contains a descriptor related to the electronic properties 
of the NCI's, Ehomo (Energy of the Highest Occupied Molecular Orbital), 
TPSA (Total Polar Surface Area) or a number of hydrogen bond acceptors. 
These models are consistent with the fact that NCI's bind at the internal 
surface of the nAChR ion channel, which is highly polar and negatively 
charged. Three of the four equations also contain a shape descriptor 



WO 2004/091492 



PCT/US2004/010978 



[- 43 -] 

(Shadow-YZ), which is consistent with the fact that the NCI's bind within 
a defined space on the receptor. Thus, the QSAR analyses describe a 
chromatographic and, as discussed above, NCI-receptor process where 
the primary driving force is electrostatic interactions between positively 
5 charged ligands and a negatively charged nAChR, which take place in the 
structurally defined central pore of the receptor. 



Table 6. QSAR equation describing affinity chromatography parameters. 



logk'= 5.255(±0.942) + 0.491(+0.092)£ HOMO + 0.01 1 8(±0.0049)7Z 
r = 0.894, s = 0.168, F = 27.929, n = 17 


Eqn. 1 


log £ on = 7.693(+0.1 1 1) - 0.00787(+0.00257)ZZ + 0.0700(±0.0237)fl&o«rf acceptors - 
-0.00276(±0.00118)rP&4 r = 0.762, s = 0.0883, F = 5.997, n = 17 


Eqn. 2 

Outlier: 

mecamylamirw 


logAr off = -3.096(±0.926)-0.454(±0.090)£ HOMO - 0.01 28(±0.00471) YZ 
r = 0.891, s = 0.165, F = 26.961, n = 17 


Eqn. 3 


logisT = 1 1 .412(±0.604) + 0.492(+0.0669)£' HOMO 
r = 0.885, s = 0.135, F = 54.130, n = 17 


Eqn. 4 







10 A 3-dimensional scatterplot of the variables associated with Eqn. 1, 

i.e. log /c 1 ., ^homo an d YZ, suggested that the whole cassette of tested NCFs 
could be subdivided into three separate clusters, Figure 13. 

The parameter k' is derived from chromatographic experiments 
using the non-linear chromatography approach described by Jozwiak 

15 et al. (Jozwiak K, Haginaka J, Moaddel R, Wainer IW. Displacement and 
nonlinear chromatographic techniques in the investigation of interaction 
of noncompetitive inhibitors with an immobilized alpha3beta4 nicotinic 
acetylcholine receptor liquid chromatographic stationary phase. Anal. 
Chem. 2002 Sep 15;74(18):4618-24) and applied in Example 2 above. 

20 Ehomo is given in electron volts (eV) and is the highest occupied molecular 
orbital energy. Ehomo is an electronic descriptor of the molecule obtained 
in molecular simulation. In the present Example Ehomo was calculated 
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using the MOPAC module in Cerius 2 software (Cerius2 v. 4.8. Accelrys 
Inc., San Diego, CA). Additional information about the Ehomo parameter 
can be found in J.M. Goodman, Chemical application of Molecular 
Modeling , c. Royal Society of Chemistry. 1998. p. 139. 
5 "Shadow YZ" is a surface area projection descriptor - the molecular 

surface is projected the YZ plane (determined by principal axis of inertia 
of the molecule) and the shadow is calculated in A 2 . In the present 
example, we used the QSAR+ module of Cerius 2 software (Cerius2 v. 4.8. 
Accelrys Inc., San Diego, CA). More information about the surface area 

10 projection descriptor can be found in Rohrbaugh et al. (Rohrbaugh RH, 
Jurs PC. Molecular shape and the prediction of high-performance liquid 
chromatographic retention indexes of polycyclic aromatic hydrocarbons. 
AnalChem. 1987 Apr 1;59(7): 1048-54). 

A cluster analysis based on the three properties was carried out 

15 using K-mean clustering method of variables and the results confirm that 
there are 3 distinct clusters. K-mean clustering is a standard clustering 
method that determines a user-specified number of clusters with the goal 
- of minimizing within-cluster- variability while maximizing between-clustery 
variability. In the present example, the method was implemented as in 

20 Statistica (STATISTICA v. 6.0. Statsoft Inc., Tulsa, OK). 

Cluster 1 was formed by four compounds (dilthiazem and 
methadone verapamil and nor-verapamil) and can be characterized by 
high values of logk' and E UOMO and YZ parameters (mean values / range: 

1.645 / 1.3 to 2.2 ; -8.93 / -9.2 to -8.6 and 64.5 / 60 to 70, respectively); 
25 Cluster 2 included 8 compounds (dextromethorphan analogs, clozapine, 
laudanosine and phencyclidine) with high values of log A;' and E UOMO but 

moderate YZ (mean values / range: 1.61 / 1.3 to 2.2; -8.64 / -9.0 to -7.7; 
and 50.0 / 45 to 60, respectively); and Cluster 3 contained 7 compounds 
(MK-801, adamantadine, bupropion, ketamine, mecamylamine, 
30 memantine, methamphetamine) with low values of log£ f , £ HOMO and YZ 
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parameters (mean value / range: 1.06 / 0.9 to 1.3; -9.45 / -9.8 to -9.1 
and 37.9 / 25 to 45, respectively). 

The analysis segregates the compounds by size and charge, with the 
smaller, more electronegative compounds appearing in Cluster 3. This 
division reflects a pharmacological reality since compounds contained in 
Cluster 3 can rapidly and deeply penetrate the luminal pore of the nAChR 
producing a high percentage of blockade per concentration of molecules. 
This would be reflected in lower IC 50 values. 

The IC 5 o values have been established for 4 in Cluster 1, 4 of the 8 
in Cluster 2 and 6 of the 8 in Cluster 3 using the Rb + efflux assay 
described above using cell lines expressing the relevant receptor. 
KXa3b4R2 is a line of human embryonic kidney 293 cells stably 
transfected with rat neuronal nicotinic acetylcholine receptor (nAChR) oc3 
and (34 subunit genes. This cell line can be obtained from Dr. Kenneth 
Kellar - Department of Pharmacology, Georgetown University, 
Washington, DC. K177 is a line of human embryonic kidney 293 cells 
stably transfected with human neuronal nicotinic acetylcholine receptor 
(nAChR) xx3 and 04 subunit genes. These cells can be obtained from Dr. " 
Daniel Bertrand, Dept. of Physiology, University of Geneva, Switzerland. 
SH-SY5Y cells are a human neuroblastoma clonal subline of the 
neuroepithelioma cell line SK-N-SH from the bone marrow. This cell ine 
can be obtained from the European Collection of Cell Cultures (ECACC), 
catalogue no. 94030304. PC- 12 Rat adrenal gland pheochromocytoma 
cells are available from the American Type Culture Collection, ATCC 
Number CRL- 1721. Results are shown in (Table 7). 

When these values were considered in relationship to the 
compounds in Clusters 2 and 3, 3 of the 4 compounds in Cluster 2 had 
IC50 values s 10 jiM while 5 of the 6 compounds in Cluster 3 had IC 50 
values < 10 juM. 
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Table 7. The IC50 values of Rb + efflux of various compounds. 



Ligand 


IC50 


Cluster 
iNumuer 


Cell line 


II1C LllcLLlVjllC 


1 . ^JLW . ^ j 


1 
± 


JV/S.IXO ptKi 


V Clcipcilllll 




1 

X 




nor- verapamil 


2.6(±1.0) 


1 


KXa3(34R2 


Dilthiazem 


2.2o(±1.0) 


1 


T/"V^, OD /I "DO 

KXaop4K2 


dextromethorphan 


0.9 (±1.1) 
1 U . 1 (± 1 . 1 U J 


0 


TV""V™ QD/1DO 
JvAXtO p 


levometnorpnan 


lU.y (±l.Uo) 


0 


TV"V™ QftADO 


dextrorphan 


^y.o (±0. / ) 




iVTvOCO p JX^i 


phencyclidine 


7.0 (+1.3) 
5.9 


0 
2 


KXa3p4R2 
SH-SY5Y 


JMLJtv-oU 1 


^20.0 (±y.o) 




iVA-CCO p -H- IS. A 


mecamylamine 


1.0 (+0.04) 


3 


KXa3p4R2 


memantine 


6.60(±0.92) 


3 


K177 (a4p2) 


amantadine 


3.44(+0.67) 


3 


K177 (a4p2) 


Bupropion 


1.4 


3 


SH-SY5Y 


Ketamine 


5.2 (±0.5) 
1.4 


3 


PC- 12 
SH-SY5Y 



The method of NCI clustering using Equation 1, above, identifies 
5 potent NCIs, i.e. those with low IC50. Compounds belonging to cluster 3 
are considered as potential NCIs and are expected to be effective in 
functional tests. Compounds in the cluster 2 are expected to express 
weaker inhibition properties. The compounds of cluster 1, which consists 
of large, bulky compounds with strong chromatographic affinity, are also 
10 expected to be potent NCIs. Initially, the IC50 value of only one of the four 
compounds in cluster 1 was known, methadone. The cluster analysis 
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predicted that verapamil, nor-verapamil and diltiazem should be effective 
NCIs of the a3p4 nAChR and functional studies confirmed this prediction. 
Functional studies were carried out using a nicotine- stimulated 
86 Rb + efflux assay on KXoc3p4R2 cell line expressing a3p4 subtype of 

5 neuronal nAChR. The studies revealed that the IC50 values of dilthiazem, 
verapamil and nor-verapamil are 2.3 |iM, 8.2|liM and 2.1 jjM respectively. 
Thus all compounds in cluster 1 are strong inhibitors. The cluster 
analysis technique is applied in Example 6 below to identify compounds 
with high potency as NCIs, which was identification was further verified 

LO by functional studies of Rb + efflux. 

The technique of cluster analysis using Eqn. 1 also suggests that 
high NCI potency can be attributed to two structurally different groups of 
compounds. It can be speculated that the two groups of compounds may 
express their inhibitory properties by two different molecular 

L 5 mechanisms . 

Example 4: Investigation of a3B2 nAChR subtype 

The 6c3134 subtype of the nAChR is "extensively characterized, easily 
accessible in stably transfected cell lines (e.g., KXa3p4R2) and widely 

20 tested in functional studies. Moreover, functional studies of this subtype 
are relatively easy. However, the transmembrane domain of cc3p4-subtype 
has some unique features not found in other subtypes. As previously 
stated, the general structure of the luminal domain is believed to be fairly 
well conserved among the subtypes of the nAChRs. However, the M2 

25 transmembrane part of p4 subunit has one critical mutation 
(phenylalanine (F) in the (34 subunit at position 15 while most other 
subunits have a valine (V) in this position (Table 1). As a result the 
nAChRs containing the p4 subunit may display significantly different 
properties than would other subunit types and may exhibit differences in 

30 the interaction of the nAChR channel with NCIs. The introduction of the 
phenylalanine moieties on the p4 subunits produces small clefts in the 
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surface of the luminal domain of the channel. These clefts play an 
important role in in the binding of NCIs, as described further below. The 
cleft explains the observed enantio selectivity between dextroemthorphan 
and levomethorphan (Jozwiak, K.; Hernandez, S.C.; Kellar, K.J.; Wainer, 
I.W. The enantioselective interactions of dextromethorphan and 
levomethorphan with the a3p4-nicotinic acetylcholine receptor: 
comparison of chromatographic and functional data. J. Chromatogr. B. 
2003, 797,423-431). 

The clefts are associated with the presence of phenylalanine in p4 
M2 domain and will not exist in other, non-p4 subtypes of the nAChR. 
Interestingly, the results from a chromatographic study which utilize an 
immobilized a3p2-nAChR column showed enantio selectivity for 
dextromethorphan and levomethorphan significantly diminished as 
compared to the immobilized a3p4-nAChR column (Table 7), further 
supporting the conclusion that the cleft is a feature of £4 subtype 
receptors that can be important for NCI activity. 



-Table 7. Comparison of enantioselectivity of DM/LM pair of enantiomers 
on two different nAChR systems. Experimental data from affinity 
chromatography - selectivity factor (a) 





oc3p4 


oc3p2 


K LM 


1.62 


1.03 


AAG° =-RT\na 


-0.29 kcal/mol 


-0.02 kcal/mol 



Based on these observations, a molecular model of oc3p2 luminal 
domain was constructed. The main difference in the structure of the p2- 
type channel is the exchange of phenylalanines from p4 helices for valines 
associated with p2 helices (See Table 1). A graphic representation of the 
model is presented in Figures 3a and 3b. Figure 3b. shows the residues 
forming the surface of the channel. The distribution of the particular rings 
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along the channel can be easily noticed. The rings are distributed as 
follows (from top to bottom): extracellular ring, leucine ring, valine ring, 
leucine ring, serine ring, threonine ring and intermediate ring (consisting 
of glutamic acid residues) and this is consistent with general 
5 considerations. The a3p2 model revealed some important differences 
when compared with the a3p4 channel. The most important is the lack of 
clefts formed on the apolar surface of the lumen. The a3(32 model is 
considered the more general of the two and represents the shape of the 
channel associated with majority of subtypes of neuronal nAChR, since 

10 there is no substantial difference in the sequence of the exposed residues 
along the channel when compared with other subunits. Only the p4 
subunit possesses a significant mutation of Val -> Phe in the valine ring. 
Therefore, the new model of the <x3p2 subtype is more homologous to 
other important subtypes of nAChR than the a3p4 model and should be 

15 considered as a general template for detailed studies of other nAChRs and 
in some perspective other members of the ligand gated ion channel 
superfamily. 

The different structure of the luminal domain of the a3p2 channel 
produces a profound change in the docking interaction of NCIs. Since 

20 there is no apolar cleft on the surface of the channel, the NCI molecules 
must find alternative interaction in the binding site. The case of special 
interest are docking simulations: dextromethorphan and levomethorphan. 
Figures 12a and 12b present the overlaid lowest energy conformations of 
these two enantiomers docked onto the model of oc3p2 luminal domain. 

25 Two molecules adopt different orientations compared to docking onto 
a3p4 model: the ligands binds primarily to the apolar part of the lumen 
with nitrogen atom interacting with serine ring. However, in contrast to 
cc3p4 docking, there is no defined cavity on the surface and only side 
interactions are possible which results in the significantly weaker energy 

50 of interaction and, what is even more important, the AG difference 
between DM-a3p2-nAChR complex and LM-a3p2-nAChR complex is 
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significantly diminished when compared to simulations on a3p4-nAChR 
model (Table 8). 

The presence of the hydrophobic cleft in the a3(34 receptor subtype 
and its absence from the a3p2 subtype presents a target for designing of 
5 compounds that are specific for one subtype over the other. 



Table 8. Enantiospecificity of dextromethorphan and levomethorphan in 
docking simulation studies of the a3p4 and a3p2 luminal channels. 





oc3p4 


a3p2 


AGdm [kcal/mol] 


-8.73 


-7.10 


AQlm [kcal/mol] 


-8.40 


-6.93 


AAG [kcal/mol] 


-0.33 


-0.17 


Example 5: Designing; c 


sf new MCI molecules 





The molecular models, clustering analysis and dynamic 



chromatographic method of the invention can be used to design molecules 
that possess enhanced activity as NCIs of nAChRs. The molecular model 
of the NCI binding site and docking studies provide an understanding of 

L5 the mechanism of non-competitive inhibition. Using the docking 
orientations of molecular NCI-nAChR complexes, we have designed 
modifications of known molecules to more strongly accommodate the 
active site and as a result obtained new compounds that express stronger 
NCI activity. Such new molecules are of interest in the pharmaceutical 

20 industry as new treatments of disorders associated with nAChR 
overactivity, e.g., as aids in smoking cessation. 

In general, the docking orientation of a putative NCI is such that 
the molecule occupies a position within the luminal channel and exhibits 
a AG of about -8.5 kcal/mol. The molecule will generally be designed to 

25 have molecular contacts with at least one, preferably 2, 3 or 4 of the side 
chains of the amino acids lining the luminal channel. Molecular contacts 



WO 2004/091492 



PCT/US2004/010978 



[" 51 -] 

that are useful in providing high binding energies (i.e. negative AG), 
include hydrogen bonds and pi orbital overlaps. 

A structure-activity relation for a NCI of a LGIC has been derived 
using the above-described methods. Thus, a compound having a bulky 
hydrophobic moiety (e.g., a phenyl or napthyl ring system or other fused 
aromatic ring system, cyclopentyl or cyclohexyl ring system, a fused ring 
system including but not limited to bicyclo [2.2.1] heptane, bicyclo [2.2.2] 
octane, morphinan and dibenzo [1.4] diazepine) and a primary, secondary 
or tertiary amino group in proximity to (i.e, approximately 5 to 10 A from, 
preferably from 5 to 8 A from, more preferably less than 7 A from) said 
hydrophobic moiety. The amino group can be directly bonded to the 
bulky hydrophobic moiety or can be linked by a spacer moiety, such as, 
but not limited to, a short hydrocarbon chain. The amino group can be 
substituted (-NR1R2, where Ri and R 2 are the same or different and are 
selected from the group consisting of H, C1-C3 alkyl, C1-C4 alkoxy, dialkyl 
keto) . The substituent is preferably one that retains a hydrogen-bonding 
potential; a preferred substituent is a keto- group, for example a dialkyl 
keto group, especially C.H 2 (C=0)CH 3 . Another preferred substituent -is- a - 
hydroxyl or alkoxyl (-CH2OH) group, e.g. ad - C 4 normal or branched 
alkoxyl group. Preferred substituted amino groups are a dialkyl keto 
amino group (e.g., HNCH 2 (C=0)CH 3 ), a hydroxyl amino group or a 
methoxy amino group. An example of such a compound is 3-methoxy- 
17-propane-2-one 9 a, 13a, 14a morphinan. 

A preferred compound designed using the method of the above 
considerations is one comprising a hydrophobic group. A preferred 
hydrophobic group comprises at least one ring that includes at least two 
conjugated unsaturated bonds, said ring optionally being fused to 
additional rings to form a ring system and said additional rings optionally 
including one or more hetero atoms. Alternatively, the hydrophobic group 
can be a hydrocarbon chain or saturated cyclic compound. The 
hydrocarbon chain can be linear or branched and preferably contains 
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from 4 to 10 carbon atoms, more preferably from 4 to 7 carbon atoms. 
The hydrocarbon chain can further include alkenyl or alkynl 
unsaturations at one or two positions. 

The compound will also preferably contain a hydrogen bond 
accepting group, which more preferably is selected from the group 
consisting of a keto group, a nitrogen-containing heterocyclic group and a 
guanidinium group. Typically, the ring or ring system and said hydrogen 
bond accepting group are joined by a linker comprising 1 to 4 carbon 
atoms and optionally containing an oxygen or sulfur atom. 

One consideration for design of effective NCI molecules is that the 
molecule will preferably span portion of the luminal domain from the 
hydrophobic region defined by the leucine and/ or valine rings to the more 
polar region defined by the serine and/ or threonine rings (Figure 2). 

The compound will preferably have activity as a non-competitive 
inhibitor of Rb + efflux of a ligand-gated neurotransmitter ion channel 
receptor with an IC50 of less than 10 juM. 

As an example of such an approach, we undertook modifications of 
dextromethorphan (DM). DM -possesses strong activity as a NCIs (ICso^lO 
M<M on a3p4-nAChR with a greatly prolonged duration of the NCI action), 
it easily passes blood-brain barrier and unlike as its enantiomer 
levomethorphan expresses little action on opioid receptors. The docking 
orientation of the DM~oc3p4~nAChR complex is presented in Figure 15. 
The molecule occupies the valine/ phenylalanine cleft on the border 
between a3-M2 helix and p4-M2 helix with an amino- group exposed for 
interaction with polar residues below this cleft. The interaction would be 
even stronger if the amino group could interact with the serine residues 
forming the serine ring, but this interaction is prevented by the distance 
of ca. 5 A separating two moieties. Study of the complex showed that the 
energy of interaction should be significantly enhanced if a methyl group 
attached to nitrogen would be exchanged into a longer moiety with 
hydrogen bond acceptors in order to allow forming strong hydrogen bonds 
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and therefore stabilizing the complex. Several possible patterns of 
dextromethorphan modification were designed and these are shown in 
Table 9. Those molecules are based on the interaction with cc3|34-nAChR 
and may be possibly selective blockers for this subtype. Compound DM- 
01 was synthesized and tested for activity in a Rb + efflux assay. The 
synthesis of DM-01 is described in Figure 14 and the data are presented 
in Figure 15. 
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Table 9. Molecules designed to inhibit a3f54-nAChR followed by AG values 
obtained in docking simulations (reference AG of dextromethorphan = - 
8.73 kcal/mol). 



Compound 


Formula 


AG 


dextrometorphan 


— o 




-8.73 kcal/mol 


DM-01 


— o 


' — — 


-9.09 kcal/mol 


DM-02 


i 

— o 




-9.35 kcal/mol 


DM-03 


I 


=o 


-10.31 kcal/mol 




\ 






DM-04 


/r-L h, t~ 
— o 




-9.39 kcal/mol 


DM-05 


<~ 

— o 


^NH 


-10.18 kcal/mol 
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Example 6: Prediction of Side Effects of Compounds Mediated by 
Non-Competitive Inhibition of Ligand-Gated Ion Channels 

The analytic methods of the present invention can also be applied to 
assessment of NCI activity of compounds compounds, both known drugs 
5 and novel compounds, to predict side effects. For example, the drugs 
verapamil, nor-verapamil and dilthiazem are commonly administered for 
treatment of high blood pressure. An undesirable side effect of these 
drugs is constipation. 

As explained above, the a3p4 nAChR subtype plays a role in 

10 regulation of gut motility and the side effects of verapamil, nor-verapamil 
and dilthiazem on gut motility have been related to NCI activity of these 
compounds against the nAChR. As an example of application of the 
analytic methods of the present invention to the investigation of drug side 
effects, we applied the cluster analysis method to predict the NCI activity 

15 of various compounds used as calcium channel blockers for treatment of 
high blood pressure, or their metabolites (MA - M6 and D-620) and listed 
in Table 10. NCI activity is predicted if the compound falls into Cluster 1. 
Predicted- NCI activity (or lack thereof) was then confirmed using the "Rh + 
efflux assay. Results are shown in Table 10. 

20 The ranges defining clusters are as above and are: Cluster 1 (low 

ICso values) log A: 1 from 1.3 to 2.2 and E UOMO from -9.2 to -8.6 and YZ from 
60 to 70; Cluster 2 (high IC 5 o values) log&' from 1.3 to 2.2 and E HOMO 
from -9.0 to -7.7 and YZ from 60 to 45 ; and Cluster 3 (low IC 5 o values) 
log A: 1 from 0.9 to 1.3; E HOMO from -9.8 to -9.1 and YZ from 45 to 25. 

25 The chromatographic method using a3p4 nAChR column was used 

to obtain experimental affinity for 13 structures (diltiazem and 5 of its 
metabolites; verapamil and 3 of its metabolites; nicardapine, nifedipine 
and amlodipine). The computational method of the invention was use to 
calculate Ehomo and YZ descriptors and the data are presented in Table 

30 10. Based on this data all compounds were assigned to respective clusters 
(Table 10). 
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After the prediction has been done based on clustering the actual 
values of IC50 were determined using the nicotine- stimulated 8 6Rb + efflux 
assay on cell line KXa3|34R2. These data are also presented in Table 10 
and the comparison of cluster method prediction with actual activity gives 
5 very good agreement - all tested ligands falls into proper categories. 

As it can be seen from Table 10 all of the tested compounds could 
be assigned to either cluster 1 or to cluster 2, which segregate them into 
two groups: very effective NCIs (IC50 < 10 jaM - cluster No. 1) and less 
effective NCIs (IC50 > 10 - cluster No. 2). The cluster analysis 
10 properly predicted the NCI activity of all 13 drugs and metabolites. 
Furthermore, the results suggest that the cardiovascular benefit 
attributed to calcium channel blocking activity may derive at least in part 
from previously unrecognized activity of inhibition of ligand-gated ion 
channels. 

15 

Table 10. Results of the cluster analysis characterization of tested calcium 
channel blockers. 



Ligand 


iog/c' 


Ehomo 


YZ 


Cluster No. 


ic 50 tmM] 


Dilthiazem 


1.64 


-8.66306 


62.29732 


1 


2.2 


MA 


1.61 


-8.6465 


66.535 


1 


4.2 


M1 


1.6 


-8.5788 


62.40288 


2 


30.4 


M2 


1.61 


-8.58 


57.4339 


2 


77.6 


M4 


1.45 


-8.4067 


61 .84228 


2 


73.2 


M6 


1.48 


-8.6465 


58.1142 


2 


63.1 


verapamil 


1.99 


-9.05746 


64.80137 


1 


8.1 


Nor-verapamil 


1.99 


-9.12446 


64.82286 


1 


2.6 


galapamil 


1.88 


-9.04879 


66.25418 


1 


6.4 


D-620 


1.25 


-9.34918 


48.52653 


2 


48.9 


nicardapine 


2.33 


-8.8397 


65.20084 


1 


2.5 


amlodipine 


2 


-8.7228 


62.93348 


1 


5.8 


nifedipine 


1.27 


-8.6323 


58.42415 


2 


24.7 



Compounds named in Table 3, Table 7 or Table 10, or specifically 
named in Figure 9 or Figure 13, and bupropion, ketamine, laudanosine, 
mecamylamine, methadone, MK-801, phenylcylclidine, ethidium, and 
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dextromethorphan are compounds known in the prior art and so are not 
considered to be inventive compounds per se within the scope of the 
present invention. Methods of the invention for non-competitively 
inhibiting a LGIC, especially a nicotinic AChR, or for treatment of a 
disease mediated by overactivity of a nicotinic AChR, exclude the use of 
bupropion, ketamine, laud ano sine, mecamylamine, methadone, MK-801, 
phenylcylclidine, ethidium, and dextromethorphan. 

The invention may be embodied in other specific forms without 
departing from the spirit or essential characteristics thereof. The foregoing 
embodiments are therefore to be considered in all respects illustrative 
rather than limiting on the invention described herein. Scope of the 
invention is thus indicated by the appended claims rather than by the 
foregoing description, and all changes that come within the meaning and 
range of equivalency of the claims are intended to be embraced therein. 

All patent and literature references cited herein are hereby 
incorporated by reference in their entirety and for all purposes, including 
the following references: 

1. ^Wainer IW, Zhang - Y, -Xiao -Yy Kellar KJ (1999) Liquid 

chromatographic studies with immobilized neuronal nicotinic 
acetylcholine receptor stationary phases: effects of receptor 
subtypes, pH and ionic strength on drug-receptor interactions. J 
Chromatogr B Biomed SciAppl 724:65-72. 

2. Zhang Y, Xiao Y, Kellar KJ, Wainer IW (1998) Immobilized nicotinic 
receptor stationary phase for on-line liquid chromatographic 
determination of drug-receptor affinities. Anal Biochem 264:22-5. 

3. Barrantes FJ. (2002) Lipid matters: nicotinic acetylcholine receptor- 
lipid interactions (Review). Mol Membr Biol 19:277-84. 

4. Morris GM, Goodsell DS, Halliday RS, et al. (1998) Automated 
docking using a Lamarckian genetic algorithm and empirical 
binding free energy function. 19:1639-62. 
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Appendix 1 AMBER Scripts for stepwise refining the model 

All Runs were made in AMBER 6.0. The computer used was SGI Octane. 
SGI Octane information is given below: 



1 195 MHZ IP30 Processor 

CPU: MIPS R 10000 Processor Chip Revision: 2.7 
FPU: MIPS R10010 Floating Point Chip Revision: 0.0 
Main memory size: 1536 Mbytes 
Xbow ASIC: Revision 1.3 
Instruction cache size: 32 Kbytes 
Data cache size: 32 Kbytes 

Secondary unified instruction/ data cache size: 1 Mbyte 
Integral SCSI controller 0: Version QL1040B (rev. 2), single ended 
Disk drive: unit 1 on SCSI controller 0 
Disk drive: unit 2 on SCSI controller 0 

Integral SCSI controller 1: Version QL1040B (rev. 2), single ended 
IOC3 serial port: ttyl 
IOC3 serial port: tty2 

IOC3 parallel port: plpl _ ^ 

Graphics board: SI 

Integral Fast Ethernet: efO, version 1, pci 2 

Iris Audio Processor: version RAD revision 12.0, number 1 



The Amber runs were made on a potassium channel receptor model that was built using the 
template structure of PDB entry 1EQ8 on Sybyl 6.8. Amber 6.0 was used to refine the structure 
that was built in sybyl 6.8. Scripts used to do the energy minimization are attached below: 



SCRIPT 1 



# 

#Script for running Sander_Classic in AMBER 6.0 

# Relaxing only Hydrogen atoms Ravi (May 23, 2002) 

# 
# 

#ALL Hs are relaxed IBELLY OPTION, \epsilon(r) 
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# 

&cntrl 

timlim=36000., imin=l, nmropt=0, 

5 ntx=l, irest=0, ntrx=l, 

ntxo=l, ntpr=10 > ntwr=0, ntwx-50, ntwv=0, ntwe=50, ntwxm=0, ntwvm=0, ntwem=0, 
ioutfm=0, ntwprt=0, 

LO ntf=l, ntb=0, idiel=0, dielc=4.0, cut=9.0, ntnb=l, nsnb=25, 
ntid~0, scnb=2.0,scee=1.2, cut2nd=0.0 

ichdna=0, 
isftrp-0, rwell=0.0, 
.5 ipol=0, 

ibelly=l, ntr=0, 

maxcyc=5000, ncyc=550, ntmin=l, dx0=0.01 ? dxm=0„05, 

0 &end 

GROUP NUMBER 1 

FIND 

* * * 

.* HI * * 
5 * HC * * 

* HP * * 

* HO * * 

* HA* * 

* HS * * 
D SEARCH 

RES 1 125 

END 

END 

SCRIPT 2 

# 

# Script for running Sander_Classic in AMBER 6.0 

► #Relaxing Hydrogen + Side-Chains Ravi (May 23, 2002) 

# 
# 
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#■ 



# Channel ALL H+SC are moving 
# 

5 &cntrl 

timlim=36000., imin-1, nmropt-0, 



ntx=l, irest=0, ntrx=l, 

10 ntxo=l, ntpr=5, ntwr=0, ntwx=50, ntwv=0, ntwe-50, ntwxm=0, ntwvm-O, ntwem=0, 
ioutfm=0, ntwprt=0, 



ntf-1, ntb-0, idiel=0, dielc=4.0, cut=9.0, ntnb=l, nsnb=25, 
ntid=0, scnb=2.0,scee=1.2, cut2nd=0.0 

L5 

ichdna=0, 
isftrp=0, rwell=0.0, 
ipol=0, 

ibelly=l, ntr=0, 

!0 

rnaxcyc=5000, ncyc=250 ? ntmin=l, dx0=0.01 ? dxm=0.05, 
&end 

.GROURN-UMBER 1 „ 
5 FIND 

* CT 3 * 

* CAB * 

* CA S * 

* OH S * 
0 * SH S * 

*S S* 
*C B* 

* N3 3 * 

* * E * 

5 SEARCH 
RES 1 125 
END 
END 

) SCRIPT 3: 



WO 2004/091492 



[- 61 -] 



PCT/US2004/0 10978 



#Script for running Sander_Classic in AMBER 6.0 
#Relaxing everything except alpha- Carbons 
#- — Ravi (May 23, 2002) 
# 

5 # 

# _ 

* Except Alpha C, all other atoms move 
# 

&cntrl 

10 timlim=36000., imin=l, nmropt=0, 

ntx=l, irest=0 } ntrx=l, 

ntxo=l, ntpr=5, ntwr=0, ntwx=50, ntwv=0, ntwe=50, ntwxm=0, ntwvm=0, ntwem=0, 
15 ioutfm=0, ntwprt=0, 

ntf=l, ntb=0, idiel=0, dielc=4.0, cut=9.0, ntnb=l, nsnb=25, 
ntid=0, scnb=2.0,scee=1.2, cut2nd=0.0 

20 ichdna=0, 

isftrp=0, rwell=0.0, 
ipol=0, 

ibelly=l, ntr=0 ? 

25 maxcyc =500 0 . ncy c=2 5 0 , n tmin= 1 , dx0=0.01, dxm=0.05 ? 
&end 

GROUP NUMBER 1 
FIND 

30 * * 3 * 

* * g * 

* * g * 

* * E * 
N N M * 

35 ccm* 

CH3 CT M * 
*HCM* 
SEARCH 
RES 1 125 
10 END 
END 
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SCRIPT 4: 

# 

5 #Script for running Sander_Classic in AMBER 6.0 

#Restrained minimization of the alpha-Carbons of the channel 

# Ravi (May 23, 2002) 

# 
# 

10 # 

# Restrained minimization of the alpha-Carbons 
# 

&cntrl 

timlim=36000. 7 imin=l, nmropt=0, 

15 

ntx=l, irest^O, ntrx= 1 , 

ntxo=l, ntpr=5 ? ntwr=0 ? ntwx=50, ntwv=0, ntwe=50, ntwxm=0, ntwvm=0, ntwem=0, 
ioutfm=0 7 ntwprt^O, 

20 

ntf=l 5 ntb«0, idiel=0, dielc=4.0, cut=9.0, ntnb=l, nsnb-25, 
ntid=0, scnb=2.0,scee=1.2, cut2nd=0.0 

ichdna=0; 
25 isftrp=0, rwell=0.0, 
ipol^O, 

ibelly=0, ntr=l, 

maxcyc=2000, ncyc=250, ntmin=l, dx0=0.01, dxm=0.05, 

30 

&end 

GROUP NUMBER 1 

10.0 

FIND 

55 ca * * * 

SEARCH 
RES 1 125 
END 
END 

[0 



#++++++++++++++++++++++++++++++ +++ + +++++ + 
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A ppendix 2 AutoGrid Parameter File 

receptor MS.pdbqs # macromolecule model 

gridfld M3.maps.fld # grid_data_file 

npts 6060120 # num. grid points in xyz 

spacing 0.375 # spacing(A) 

gridcenter 0.009 0.026 -0.172 # xyz- coordinates or auto 

types CANOH # atom type names 

smooth 0.5 # store minimum energy w/in rad(A) 

map M3.C.map # atom-specific affinity map 



nbp_r_eps 


4.00 0.0222750 


12 


6 


# C-C lj 


nbp_r_eps 


3.75 0.0230026 


12 


6 


# C-N lj 


nbp_r_eps 


3.60 0.0257202 


12 


6 


# C-O lj 


nbp_r_eps 


4.00 0.0257202 


12 


6 


# C-S lj 


nbp_r_eps 


3.00 0.0081378 


12 


6 


# C-H lj 


nbp„r_eps 


3.00 0.0081378 


12 


6 


# C-H lj 


nbp_r_eps 


3.00 0.0081378 


12 


6 


# C-H lj 



sol_par 12.77 0.6844 # C atomic fragmental volume, solvation parameters 

constant 0.000 # C grid map constant energy 

map M3.A.map # atom-specific affinity map 



nbp_ 


jr_ 


_eps 


4.00 0.0222750 


12 


6 


# A-C lj 


nbp_ 


jr_ 


_eps 


3.75 0.0230026 


12 


6 


# A-N lj 


nbp_ 


_r_ 


_eps 


3.60 0.0257202 


12 


6 


# A-O lj 


nbp. 


_r_ 


_eps 


4.00 0.0257202 


12 


6 


# A-S lj 


nbp_ 


X. 


_eps 


3.00 0.0081378 


12 


6 


# A-H lj 


nbp_ 


T_ 


_eps 


3.00 0.0081378 


12 


6 


# A-H lj 


nbp_ 


r_ 


_eps 


3.00 0.0081378 


12 


6 


# A-H lj 



sol_par 10.80 0.1027 # A atomic fragmental volume, solvation parameters 

constant 0.000 # A grid map constant energy 

map M3.N.map # atom-specific affinity map 



nbp_ 


_r_ 


_eps 


3.75 0.0230026 12 


6 


# N-C lj 


nbp_ 


x. 


_eps 


3.50 0.0237600 12 


6 


# N-N lj 


nbp_ 


_r_ 


_eps 


3.35 0.0265667 12 


6 


# N-O lj 


nbp_ 


_r_ 


_eps 


3.75 0.0265667 12 


6 


# N-S lj 


nbp_ 


_r_ 


_eps 


1.90 0.3280000 12 


10 


# N-H hb 


nbp_ 


_r_ 


_eps 


1.90 0.3280000 12 


10 


# N-H hb 


nbp_ 


_r_ 


_eps 


1.90 0.3280000 12 


10 


# N-H hb 



sol_par 0.00 0.0000 # N atomic fragmental volume, solvation parameters 

constant 0.000 # N grid map constant energy 

map MS.O.map # atom-specific affinity map 
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nbp__r_eps 3.60 0.0257202 12 6 # O-C lj 
nbp_r__eps 3.35 0.0265667 12 6 # O-N lj 
nbp_r_eps 3.20 0.0297000 12 6 # O-O lj 
nbp_r_eps 3.60 0.0297000 12 6 # OS lj 
nbp_r_eps 1.90 0.3280000 12 10 # O-H hb 
nbp_r_eps 1.90 0.3280000 12 10 # O-H hb 
nbp_r.eps 1.90 0.3280000 12 10 # O-H hb 

sol_par 0.00 0.0000 # O atomic fragmental volume, solvation parameters 

constant 0.236 # O grid map constant energy 

map M3.H.map # atom-specific affinity map 

nbp__r_eps 3.00 0.0081378 12 6 # H-C lj 

nbp_r_eps 1.90 0.3280000 12 10 # H-N hb 

nbp_jr_eps 1.90 0.3280000 12 10 # H-O hb 

nbp_r_eps 3.00 0.0093852 12 6 # H-S lj 

nbp_r_eps 2.00 0.0029700 12 6 # H-H lj 

nbp_r_eps 2.00 0.0029700 12 6 # H-H lj 

nbp_r_eps 2.00 0.0029700 12 6 # H-H lj 

sol_par 0.00 0.0000 # H atomic fragmental volume, solvation parameters 

constant 0.1 18 # H grid map constant energy 

elecmap M3.e.map # electrostatic potential map 

dielectric 15.0 # <0, distance-dep.diel;>0, constant 

fmap M3.f.map # floating point potential gridmap 



Appendix 3 AutoBock Parameter File 



seed pid time 

types CANOH 

fid M3.maps.fld 

map M3.C.map 

map MS.A.map 

map MS.N.map 

map M3.0.map 

map MS.H.map 

map MS.e.map 

move DMT.out.pdbq 

about -0.088 0.126 0.069 

tranO random 

quatO random 

ndihe 1 

diheO random 

tstep 2.0 



# seeds for random generator 

# atom type names 

# grid_data_file 

# atom-specific affinity map 

# atom-specific affinity map 

# atom-specific affinity map 

# atom-specific affinity map 

# atom-specific affinity map 

# electrostatics map 
# small molecule 

# small molecule center 

# initial coordinates/ A or random 

# initial quaternion 

# number of active torsions 

# initial dihedrals (relative) or random 

# translation step/ A 
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qstep 50.0 # quaternion step/deg 

dstep 50.0 # torsion step/deg 

torsdof 10.3113 # torsional degrees of freedom and coeffiecent 

intnbp_r_eps 4.00 0.0222750 12 6 # C-C Ij 
5 intnbp_r_eps 4.00 0.0222750 12 6 # C-A lj 

intnbp_r_eps 3.75 0.0230026 12 6 # C-N lj 

intnbp_r_eps 3.60 0.0257202 12 6 # C-O lj 

intnbp_r_eps 3.00 0.0081378 12 6 # C-H lj 

intnbp„r„eps 4.00 0.0222750 12 6 # A-A lj 
10 intnbp_r_eps 3.75 0.0230026 12 6 # A-N lj 

intnbp_r_eps 3.60 0.0257202 12 6 # A-O lj 

intnbp_r_eps 3.00 0.0081378 12 6 # A-H lj 

intnbp_r_eps 3.50 0.0237600 12 6 # N-N lj 

intnbp„r_eps 3.35 0.0265667 12 6 # N-O lj 
15 intnbp_r_eps 2.75 0.0084051 12 6 # N-H lj 

intnbp_r„eps 3.20 0.0297000 12 6 # O-O lj 

intnbp_r_eps 2.60 0.0093852 12 6 # O-H lj 

intnbp_r_eps 2.00 0.0029700 12 6 # H-H lj 

outlev 1 # diagnostic output level 

20 rmstol 0.5 # cluster_tolerance /A 

extnrg 1000.0 # external grid energy 

eOmax 0.0 10000 # max initial ernergy; max number of retries 

ga_pop_size 50 # number of individuals in population 

-ga_num_evals 5000000 # maximum number of energy evaluations 

15 ga_num_generations 27000 # maximum number of generations 

ga_elitism 1 # number of top individuals to survive to next generation 

ga_mutation_rate 0.02 # rate of gene mutation 

ga_crossover_rate 0.8 # rate of crossover 

ga_window_size 10 # 
50 ga_cauchy_alpha 0.0 # Alpha parameter of Cauchy distribution 

ga_cauchy_beta 1.0 # Beta parameter Cauchy distribution 

set_ga # set the above parameters for GA or LGA 

sw_max_its 300 # iterations of Solis & Wets local search 

sw_max_succ 4 # consecutive successes before changing rho 

»5 sw_max_fail 4 # consecutive failures before changing rho 

sw_rho 1.0 # size of local search space to sample 

sw_lb_rho 0.01 # lower bound on rho 

ls__search_freq 0.06 # probability of performing local search on individual 

set_pswl # set the above pseudo-Solis & Wets parameters 

0 ga_run 50 # do this many hybrid GA-LS runs 

analysis # perform a ranked cluster analysis 
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Appendix 4 : Atomic Coordinates of the Luminal Channel of a 

cc3p4 nAChR Ion Channel 

pdb file of the <x3p4 nAChR model 
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CLAIMS 

What is claimed is: 

5 1 . A computer system comprising: 

i) a memory storing positional data of the atomic coordinates of the 
transmembrane portion of at least one subunit of a ligand-gated 
neurotransmitter receptor protein; and 

ii) a processor generating a molecular model having a three 
10 dimensional shape representative of the luminal domain portion of the 

ligand-gated neurotransmitter receptor based the positional data. 

2. The computer system of claim 1 ? wherein the memory stores data of 
the atomic coordinates of at least an a3 chain of a nicotinic acetylcholine 

15 receptor and a (34 chain of a nicotinic acetylcholine receptor. 

3. The computer system of claim 1, wherein the memory stores data of 

the atomie coordinates - of at least one polypeptide having an amino acid - — 
sequence selected from the group consisting of SEQ ID NOS: 1 through 
20 15. 

4. The computer system of claim 1 ? wherein the processor generates a 
molecular model of the luminal domain portion of a ligand-gated 
neurotransmitter receptor having a subunit stoichiometry ranging from 

25 (a) 5 (0)o to (a) 2 (p) 3 , or (a) 2 p§y. 

5. The computer system of claim 4, in which the stoichiometry is (a)2(p)3, 

6. The computer system of claim 5, in which the a subunits are a3 
30 subtype and the (3 subunits are (32 or (34 subtype. 
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7. The computer system of claim 1, wherein the data comprises the 
atomic coordinates of a portion of the transmembrane portion of the 
subunit consisting of the amino acid sequence of residues 8 to 19 of at 
least one of SEQ ID NOS: 1-15. 

5 

8. The computer system of claim 1, further comprising a removable, 
computer readable medium in which the data are stored. 

9. A method for screening a compound for activity as a non-competitive 
10 inhibitor of a ligand-gated neurotransmitter receptor comprising; 

identifying as a compound having activity as a non-competitive inhibitor 
as one having a AG less than -6 kcal/mol by docking a model of the 
compound in a model of the three-dimensional shape of the luminal 
portion of the ligand-gated neurotransmitter receptor. 

15 

10. The method of claim 9, in which the docking step is performed using 
a computer program implementing a genetic algorithm. 

11. The method of claim 9 ? further comprising assaying the binding of 
20 the compound to the receptor and confirming as a non-competitive 

inhibitor one having a k' of greater than 8, the assaying step being 
performed either before or after the docking step. 

12. The method of claim 11, further comprising modifying the structure 
25 of the compound confirmed as a non-competitive inhibitor and repeating 

the method using the compound having the modified structure. 

13. A method for screening a compound for activity as a non- 
competitive inhibitor of a ligand-gated neurotransmitter receptor 

30 comprising: 
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i) identifying as a compound having activity as a non-competitive 
inhibitor as one having a AG less than -6 kcal/mol by docking a model of 
the compound in a model of the three-dimensional shape of the luminal 
domain portion of the ligand-gated neurotransmitter receptor; 
5 ii) testing the compound for binding to the ligand-gated 

neurotransmitter receptor and/ or inhibiting the ion-channel activity 
thereof; 

iii) obtaining as a non-competitive inhibitor of the ligand-gated 
neurotransmitter receptor the compound that binds to the ligand-gated 
10 neurotransmitter receptor with a k ? value greater than 8 and being 
displaced in chromatographic experiments by mecamylamine and/ or 
inhibits the ion channel activity of the ligand-gated neurotransmitter 
receptor in nicotine stimulated 86 Rb + efflux with IC50 lower than 100 pM. 

15 14. A method for making a pharmaceutical composition comprising 
admixing a compound identified by the method of claim 13 and a 
pharmaceutical^ acceptable carrier. 

15. The method of claim 13, in which the ligand-gated neurotransmitter 
20 receptor is a neuronal receptor. 

16. The method of claim 13, in which the ligand-gated neurotransmitter 
receptor is a muscular receptor. 

25 17. A compound comprising: a hydrophobic moiety and primary, 
secondary or tertiary amino group located 5 to 10 A from said 
hydrophobic moiety, said compound having activity of inhibiting the ion- 
channel activity of a ligand-gated neurotransmitter receptor in an assay of 
nicotine stimulated 86 Rb + efflux with an IC50 lower than 50 jjM. 

30 
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18. The compound of claim 17, wherein the hydrophobic moiety is 
selected from the group consisting of an alkyl chain, a phenyl ring, a 
napthyl ring, a cyclopentyl ring, a cyclohexyl ring, and a fused ring. 

5 19. The compound of claim 18, wherein the hydrophobic moiety is a 
fused ring selected from the group consisting of bicyclo [2.2.1] heptane, 
bicyclo [2.2.2] octane, morphinan and dibenzo [1.4] diazepine. 

20. The compound of claim 17, that does not compete with the 
10 neurotransmitter ligand of a receptor for binding of a neurotransmitter 
ligand binding site of the receptor located on the external surface at the 
interface of two subunits in a pocket approximately 30-35 A from the 
transmembrane portion of the subunits. 

15 21. The compound of claim 17, wherein said receptor is a neuronal 
nicotinic acetylcholine receptor. 

22. The compound of -claim 17 , wherein the receptor is a muscular™ — 
nicotinic acetylcholine receptor. 

20 

23. A compound comprising a hydrophobic moiety and a primary, 
secondary or tertiary amino group placed 5 to 10 A from said hydrophobic 
moiety, said compound having activity of binding to the luminal domain 
site of a ligand-gated neurotransmitter receptor with a AG of at least -6 

25 kcal/mol. 

24. The compound of claim 23, that does not compete with the 
neurotransmitter ligand of the receptor for binding to the 
neurotransmitter ligand binding site of the receptor. 

30 
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25. The compound of claim 23, wherein said receptor is a neuronal 
nicotinic acetylcholine receptor. 

26. The compound of claim 23, wherein the receptor is a muscular 
5 nicotinic acetylcholine receptor. 

27. A kit comprising a computer readable medium having stored therein 
positional data of the atomic coordinates of the transmembrane portion of 
at least one subunit of a ligand-gated neurotransmitter receptor protein; 

10 and 

a composition comprising the ligand-gated neurotransmitter 
receptor protein. 

28. The kit of claim 27, in which the composition comprising the receptor 
15 protein is an affinity chromatography medium. 

29. The kit of claim 28, in which the receptor protein is a nicotinic 

-acetylcholine receptor protein. 

20 30. A computer readable medium having stored therein positional data of 
the atomic coordinates of the transmembrane portion of at least one 
subunit of a ligand-gated neurotransmitter receptor protein. 

31. The computer readable medium of claim 30, wherein the data 
25 comprise the atomic coordinates of at least an a3 chain of a nicotinic 

acetylcholine receptor and a P4 chain of a nicotinic acetylcholine receptor. 

32. The computer readable medium of claim 30, wherein the data 
comprise the atomic coordinates of at least one polypeptide having an 

30 amino acid sequence selected from the group consisting of SEQ ID NOS: 1 
through 15. 
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33. The computer readable medium of claim 30, wherein the data 
comprise a molecular model of the luminal domain portion of a ligand- 
gated neurotransmitter receptor having a subunit stoichiometry ranging 

5 from (a) 5 (p)o to (a) 2 (p)3, or (a) 2 p5y. 

34. The computer readable medium of claim 33, in which the 
stoichiometry is (0)20)3. 

10 35. A method for non~competitively inhibiting a ligand-gated ion channel 
receptor comprising contacting the ligand-gated ion channel receptor with 
a compound of claim 17, except for bupropion, ketamine, laudanosine, 
mecamylamine, methadone, MK-801, phenylcylclidine , ethidium, and 
dextromethorphan. 

15 

36. A method for treating Tourette's syndrome, schizophrenia, a 
cognitive disorder, pain, anxiety, depression, neuro degeneration or an 
addiction caused by an overactive ligand-gated—ionr channel receptor,- 
comprising administering to a subject an amount of a compound of claim 

20 16 effective to inhibit ion flux through said ligand-gated ion channel. 

37. The method of claim 36, that is a method of treating smoking 
addiction and the receptor is a neuronal nicotinic acetylcholine receptor. 

25 38. A method for evaluating cardiovascular toxicity or effects on gut 
motility of a compound comprising testing a compound for activity in the 
method of claim 13, wherein the receptor is a neuronal nicotinic 
acetylcholine receptor subtype. 

30 39. A method for identifying a compound as a non-competitive inhibitor 
of a ligand-gated neurotransmitter ion channel receptor comprising: 



WO 2004/091492 



PCT/US2004/0 10978 



[- 107 -] 

i) determining the log k' for binding of the compound to the ligand- 
gated neurotransmitter ion channel receptor, or to a subunit thereof; 

ii) determining the energy of the highest occupied molecular orbital 
of the compound in eV; 

5 iii) determining the area in A 2 of a plane projection of the 

compound; and 

iv) identifying as an effective non-competitive inhibitor a compound 
exhibiting all of log k' from 1.1 to 1.9, an energy of the highest occupied 
molecular orbital from -8.6 to -9.2 eV and a YZ shadow of 25 to 50 A 2 . 

10 

40. A method for making a pharmaceuctical composition comprising 
admixing a compound identified by the method of claim 39 with a 
pharmaceutically acceptable carrier. 

15 41. The method of claim 39, in which k' is determined using docking of a 
computational model of the compound to a computational model of the 
luminal channel of the ligand-gated neurotransmitter ion channel 
receptor.. 

20 42. The method of claim 39, in which k ? is determined by 
chromatography of the compound on an affinity matrix comprising at 
least one subunit polypeptide of the ligand-gated neurotransmitter ion 
channel receptor. 

25 43. The method of claim 39, in which the ligand-gated neurotransmitter 
receptor has a subunit stoichiometry ranging from (a)sO)o to (a)2(P)3, or 
(a)2PSy. 



44. The method of claim 43, in which the stoichiometry is (a)2(p)3, 

30 
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45. The method of claim 44, in which the model of the luminal channel 
of the ligand-gated neurotransmitter ion channel receptor has the atomic 
coordinates of the a3(34 subtype receptor or the a3(52 subtype receptor 
shown in Appendix 4 or Appendix 5, respectively. 

5 

46. A compound that is a derivative of dextromethorphan having the 

I. 

nitrogen-bound methyl group substituted by a Ci~6 alkyl group bearing a 
hydrogen-bond accepting group. 

10 47. The compound of claim 46, in which the hydrogen bond accepting 
group is a keto group, a guanidinium group or a nitrogen-containing 
heterocyclic group. 

48. The compound of claim 47, in which the nitrogen-containing 
15 heterocyclic group is a pyrollidine, imidazolidine, piperidine, 
hexahydropyrimidine or pyrimidine group. 

-49. A compound comprising 

a hydrophobic group comprising a saturated or unsaturated alkyl 
20 chain containing 4 to 10 carbon atoms, a saturated hydrocarbon ring 
containing 5 or 6 carbon atoms, or at least one ring that includes at least 
two conjugated unsaturated bonds, said ring optionally being fused to 
additional rings to form a ring system and said additional rings optionally 
including one or more hetero atoms; 
25 a hydrogen bond accepting group selected from the group 

consisting of a keto group, a nitrogen-containing heterocyclic group and a 
guanidinium group; 

a linker joining said hydrophobic group and said hydrogen bond 
accepting group and comprising 1 to 4 carbon atoms and optionally 
30 containing an oxygen or sulfur atom; 
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the compound having activity as a non-competitive inhibitor of Rb + efflux 
of a ligand-gated neurotransmitter ion channel receptor with an IC50 of 
less than 10 juM. 

50. The compound of claim 49 , in which the ring is a planar, aromatic 
ring system. 

51. The compound of claim 49, in which the ring is selected from the 
group consisting of a phenyl ring, a napthyl ring, morphinan and dibenzo 
[1.4] diazepine. 
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Figure 6 
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Figure 8 
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Figure 13 
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2004-04-09 Sequence Listing as filed 
SEQUENCE LISTING 

<110> Irving w. WAINER et al . 

<120> COMPUTER-BASED MODEL FOR IDENTIFICATION AND CHARACTERIZATION OF 
NON-COMPETITIVE INHIBITORS OF NICOTINIC ACETYLCHOLINE RECEPTORS 
AND RELATED LIGAND-GATED ION CHANNEL RECEPTORS 

<130> 1173-1025PUS2 

<160> 15 

<170> Patentln version 3.2 

<210> 1 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi sc_feature 
<222> (1) . . (23) 

<223> Table 1 Delta Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 1 

Glu Lys Met ser Thr Ala lie ser val Leu Leu Ala Gly Ala val Phe 
15 10 15 

Leu Leu Leu Thr ser Gly Arq 
20 



<210> 2 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi scfeature 
<222> (1). .(23) 

<223> Table 1 Gamma Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 2 

Gin Lys cys Thr Leu ser lie ser val Leu Leu Ala Gin Thr lie Phe 
15 10 15 

Leu Phe Leu lie Ala Gin Lys 
20 



<210> 3 

<211> 23 

<212> PRT 

<213> unsure 

<220> 

<221> mi sc_f eature 

<222> (1) . . (23) 

<223> Table 1 Alpha 1 sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 3 
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2004-04-09 Sequence Listing as filed 
Glu Lys Met Thr Leu ser lie Ser val Leu Leu ser Leu Thr val phe 
15 10 15 

Leu Leu val He val Glu Leu 
20 



<210> 4 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi s cofeature 
<222> CI) . . (23) 

<223> Table 1 Alpha 3 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 4 

Glu Lys val Thr Leu Cys lie Ser Val Leu Leu Ser Leu Thr Val Phe 
1 5 10 15 

Leu Leu val lie Thr Glu Thr 
20 



<210> 5 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi sc_f eat u re 
<222> (1) . . C23) 

<223> Table 1 Alpha 4 sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 5 

Glu Lys lie Thr Leu Cys lie Ser val Leu Leu ser Leu Thr Val Phe 
1 5 10 15 

Leu Leu Leu lie Thr Glu lie 

20 



<210> 6 

<211> 23 

<212> PRT 

<213> unsure 

<220> 

<221> misc_feature 
<222> CD - - C23) 

<223> Table 1 Alpha 5 sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 6 

Glu Lys lie Cys Leu cys Thr ser Val Leu Val ser Leu Thr val Phe 
15 10 15 

Leu Leu Val lie Glu Glu lie 
20 
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<210> 7 
<211> 23 
<212> PRT 
<213> Unsure 

<220> 

<221> mi sc_f eature 
<222> CI) - ■ (23) 

<223> Table 1 Alpha 6 sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 7 

Glu Lys val Thr Leu cys lie ser val Leu Leu Ser Leu Thr Val Phe 
1 5 10 15 

Leu Leu Val lie Thr Glu Thr 
20 



<210> 8 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi sc_f eature 
<222> (1) . . (23) 

<223> Table 1 Alpha 7 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 8 

Glu Lys lie Ser Leu Gly lie Thr val Leu Leu ser Leu Thr val Phe 
1 5 10 15 

Met Leu Leu Val Ala Glu lie 

20 



<210> 9 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi sc_f eature 
<222> (1) . . (23) 

<223> Table 1 Alpha 9 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 9 

Glu Lys val ser Leu Gly Val Thr lie Leu Leu Ala Met Thr Val Phe 
15 10 15 

Gin Leu Met val Ala Glu lie 
20 



<210> 10 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi sc_f eature 
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<222> CD - - (23) 

<223> Table 1 Alpha 10 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 10 

Glu Lys Val Ser Leu Gly Val Thr Val Leu Leu Ala Leu Thr val Phe 
15 10 15 

Gin Leu lie Leu Ala Glu ser 
20 



<210> 11 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> mi s cofeature 
<222> (1) . . C23) 

<223> Table 1 Beta 1 sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 11 

Glu Lys Met Gly Leu Ser lie Phe Ala Leu Leu Thr Leu Thr val Phe 
15 10 15 

Leu Leu Leu Leu Ala Asp Lys 
20 



<210> 12 

<211> 23 

<212> PRT 

<213> Unsure 

<220>~ 

<221> mi sc_f eat u re 
<222> CI) - - C23) 

<223> Table 1 Beta 2 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 12 

Glu Lys Met Thr Leu Cys lie Ser val Leu Leu Ala Leu Thr Val Phe 
1 5 10 15 

Leu Leu Leu lie ser Lys lie 
20 



<210> 13 

<211> 23 

<212> PRT 

<213> Unsure 

<220> 

<221> misc_feature 
<222> CD - -C23) 

<223> Table 1 Beta 3 Sequence - Transmembrane domain of ligand gated ion 
channel subunit 

<400> 13 

Glu Lys Leu ser Leu Ser Thr Ser Val Leu val ser Leu Thr val Phe 
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2004-04-09 Sequence Listing as filed 
1 5 10 15 

Leu Leu val lie Glu Glu lie 
20 

<210> 14 
<211> 23 
<212> PRT 
<213> Unsure 

<220> 

<221> misc_feature 
<222> (I) . . (23) 

<223> Table 1 Beta 4 Sequence - Transmembrane domain of liqand gated ion 
channel subunit 

<400> 14 

Glu Lys Met Thr Leu cys lie Ser val Leu Leu Ala Leu Thr Phe Phe 
1 5 10 15 

Leu Leu Leu lie Ser Lys lie 
20 

<210> 15 
<211> 23 
<212> PRT 
<213> Unsure 

<220> 

<221> mi sc_f eat u re 
<222> (1). .(23) 

<223> Table 1 Epsilon sequence - Transmembrane domain of liqand gated ion 
channel subunit 

<400> 15 

Gin Lys cys Thr~val ser lie Asn val Leu Leu Ala Gin Thr val Phe 
1 5 10 15 

Leu Phe Phe Leu lie Ala Gin 
20 

5 of 5 



Page 5 



